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ABSTRACT  OF  THE  THESIS 

Modeling  Present  and  Future  River  Runoff  Using 
Global  Atmospheric  Models 

by  SCOTT  C  VAN  BLARCUM 

Thesis  Director: 

Dr  James  Miller 

A  global  atmospheric  model  is  used  to  calculate  the  monthly  river  runoff  for  30  of 
the  world's  major  rivers  for  the  present  climate  and  for  a  doubled  CO2  climate.  The  model 
has  a  horizontal  resolution  of  4°  X  5°,  but  the  runoff  from  each  model  box  is  quartered 
and  added  to  the  appropriate  river  drainage  basin  on  a  2°  X  2.5°  resolution  A  new 
routing  scheme  is  used  to  allow  runoff  calculated  for  a  particular  grid  box  to  flow  to  an 
adjacent  downstream  grid  box  and  ultimately  to  the  mouth  of  the  river  The  total 
instantaneous  runoff  leads  runoff  at  the  mouth  by  one  to  two  months.  The  model¬ 
generated  runoff  at  the  mouth  is  compared  to  observations  for  several  different 
simulations  The  runoff  peaks  of  high-latitude  rivers  are  due  to  spring  snow  melt  and 
there  is  a  time  lag  between  when  the  snow  melts  and  when  the  melt  water  reaches  the 
mouth.  The  new  routing  scheme  allows  the  calculation  of  runoff  at  any  location  in  the 
river  basin.  Model-generated  river  runoff  and  precipitation  for  the  Mississippi  River  and 
its  tributaries  are  analyzed  for  the  present  climate,  where  annual  precipitation  is  within  5% 
of  the  observed  precipitation.  However,  model-generated  monthly  precipitation  is  too 
high  in  the  spring  and  too  low  in  the  summer  and  fall.  In  a  model  simulation  with  doubled 
CO2,  river  runoff  increases  for  27  of  the  30  rivers  and  in  most  cases  coincides  with 
increased  precipitation.  All  high-latitude  rivers  show  an  increase  in  precipitation  and 
runoff  with  a  shift  in  the  runoff  maximum,  approximately  one  month  earlier,  due  to  an 
earlier  snow  melt  season.  In  a  doubled  CO2  climate,  snow  mass  decreases  for  mid  and 
high-latitude  rivers  in  North  America  and  northwestern  Asia,  but  increases  for  rivers  in 
northeastern  Asia,  where  observed  winter  temperatures  average  -30°  to  -50°  C 
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Chapter  1.  Introduction 

Global  warming  associated  with  the  increase  in  greenhouse  gases,  such  as  carbon 
dioxide  (C02),  continues  to  be  a  major  topic  within  the  scientific  and  political  fields  Most 
people  associate  temperature  change  with  global  warming,  however,  the  earth's  hydrologic 
cycle  will  also  be  affected  by  the  increase  of  C02  Over  land,  the  major  components  of  the 
hydrologic  cycle  include  precipitation,  ground  water,  evapotranspiration,  ice,  and  runoff 
Runoff  is  generated  by  precipitation  and  snow  melt  Because  evaporation  exceeds 
precipitation  over  the  oceans,  it  is  runoff  from  the  land  to  the  ocean  that  keeps  the 
hydrologic  cycle  in  equilibrium.  River  runoff  depends  on  precipitation,  moisture  storage 
within  the  soil  (dependent  on  evaporation,  soil  type,  and  vegetation  cover),  and  flow  rate 
due  to  topography  [Thomtfmaite,  1955],  The  world’s  20  largest  rivers  account  for  about 
40%  of  the  total  continental  runoff  [Baumgartner  and  Reichel ,  1 975] 

Several  atmospheric  general  circulation  models  (GCMs)  have  been  developed  to 
allow  detailed  scientific  studies  of  the  hydrologic  cycle  [Hansen  et  a/.,  1983,  Gordon  and 
Stern ,  1982,  Manabe  and  Hahn ,  1981,  Pitcher  et  al ,  1983;  and  Williamson ,  1983] 
Russell  and  Miller  [  1 990]  compared  mean  annual  river  runoff  from  a  four-year  GCM 
simulation  with  observations  Kuhl  and  Miller  [1992]  extended  that  work  to  examine 
monthly  river  runoff  statistics  for  wet,  dry,  and  moderately  wet  river  basins  Potential 
changes  in  the  hydrologic  cycle  (specifically  river  runoff)  due  to  global  warming  will  affect 
many  areas,  including  agriculture,  water  resources,  and  hydroelectric  power 

Potential  changes  in  river  runoff  have  caused  hydrologists  to  generally  focus  on 
calculating  river  runoff  at  smaller  scales  (about  1  kilometer)  than  those  considered  by 
global  atmospheric  modelers.  Most  studies  are  centered  around  specific  river  basins  and 
assess  the  potential  changes  in  river  runoff  due  to  global  warming  Vordsmarty  et  al 
[1989]  have  developed  a  continental  drainage  basin  runoff  scheme  and  applied  it  to  the 
Amazon  River  basin  using  observed  data  as  input.  This  is  one  of  the  first  attempts  to 
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model  river  runoff  over  much  larger  scales  The  model  allows  the  calculation  of  water 
volume  from  one  grid  box  to  another,  however,  the  model  is  based  on  observed  data  and 
can  not  be  used  directly  to  calculate  potential  changes  in  nver  runoff  due  to  global 
warming 

There  are  two  primary  methods  used  to  study  effects  of  global  warming  on  river 
runoff  using  GCMs.  The  first  method  uses  temperature  and  precipitation  variables  from 
doubled  CO2  simulations  as  input  to  a  hydrologic  model  to  produce  runoff  statistics.  This 
is  referred  to  as  an  off-line  method.  Regional  river  runoff  changes  due  to  global  warming 
using  the  off-line  method  are  discussed  in  Gleick  [1987]  and  I.ettenmaier  and  Gan  [1990] 
The  off-line  method  can  also  be  used  to  input  temperature  and  precipitation  based  on 
scenarios  from  other  models  and  studies.  Several  scenario  studies  have  been  developed 
which  focus  on  specific  river  basins  around  the  world  [Ayers  el  a/.,  1990,  Flasehko  et  al., 
1987,  and  Roos ,  1989],  The  second  method,  which  computes  river  runoff  as  part  of  the 
GCM  simulation,  is  referred  to  as  an  on-line  method  Changes  in  river  runoff  using  this 
method  are  discussed  in  Miller  and  Russell  [1992],  Rind  et  al.  [1990]  analyze  both  on¬ 
line  and  off-line  methods  with  respect  to  drought  conditions  using  potential 
evapotranspiration.  They  suggest  that  drought  intensification  is  understated  in  most 
GCM  simulations  because  they  lack  realistic  land  surface  parameterizations. 

The  purpose  of  this  thesis  is  to  obtain  a  better  understanding  of  land  surface 
hydrology,  particularly  river  runoff,  as  it  is  modeled  in  global  atmospheric  GCMs  and 
examine  potential  changes  in  river  runoff  due  to  global  warming  using  the  on-line  method 
with  river  runoff  calculated  from  the  Goddard  Institute  for  Space  Studies  (GISS)  GCM  of 
Hansen  et  al.  [1983]  This  study  will  utilize  the  version  of  the  GISS  model  which  has  a 
horizontal  resolution  of  4°  latitude  by  5°  longitude  and  nine  vertical  layers.  River  runoff 
is  calculated  from  several  different  three  to  five  year  simulations.  Each  of  the  4°  X  5°  grid 
boxes  overlays  four  2°  X  2.5°  grid  boxes  which  are  assigned  to  river  basins.  The  observed 
river  runoff  is  defined  as  the  flow  rate  at  a  certain  location  within  the  river  basin.  A  new 
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river  routing  scheme  is  used  to  obtain  runoff  at  various  sites  within  the  river  basin  A 
description  of  the  hydrology  scheme  used  in  the  model  and  the  new  routing  scheme  will  be 
given  The  objectives  follow: 

1.  Apply  a  new  river  routing  scheme  with  the  GISS  two-layer  grid  point 
hydrology  scheme  and  examine  the  monthly  runoff  for  approximately  30 
major  river  basins  In  particular,  high  latitude  river  basins  will  be  examined 
to  determine  how  snow  melt  affects  monthly  runoff  and  dry  river  basins 
will  be  examined  to  leam  why  the  runoff  is  poorly  simulated  by  the  model. 

2.  The  new  routing  scheme  which  allows  flow  from  one  grid  box  to 
another  will  be  used  to  simulate  the  runoff  of  the  Mississippi  river  basin  by 
examining  the  runoff  of  its  major  tributaries  (Missouri,  Illinois,  Arkansas, 
and  Ohio  rivers). 

3.  Examine  how  monthly  river  runoff  responds  to  the  doubling  of 
atmospheric  carbon  dioxide  by  examining  model  simulations  with  an 
emphasis  on  high  latitude  rivers  dominated  by  spring  snow  melt 

River  runoff  and  precipitation  from  approximately  30  river  basins  throughout  the 
world  representing  dry,  wet,  and  moderately  wet  climates  as  well  as  high,  middle,  and  low 
latitudes  will  be  examined  based  on  several  GCM  simulations  of  the  GISS  model  Results 
discussed  in  Russell  and  Miller  [1990],  Miller  and  Russell  [1992],  and  Kuhl  and  Miller 
[1992]  will  be  extended  by  allowing  river  runoff  to  be  routed  between  grid  boxes  within  a 
river's  drainage  basin.  These  simulations  will  be  used  to  analyze  the  similarities  and 
differences  between  monthly  runoff  calculated  with  this  new  routing  scheme  and 
compared  with  the  total  monthly  runoff  into  a  river  basin  as  in  Kuhl  and  Miller  [1992], 

With  this  routing  scheme  in  place,  the  components  of  the  Mississippi  river  basin 
will  be  studied  separately  These  include  the  upper  Mississippi  River,  Illinois  River, 
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Missouri  River,  Ohio  River,  and  Arkansas  River  The  basin  will  be  separated  into  2°  X 
2.5°  grid  boxes  The  Mississippi  River  and  its  tributaries  will  be  examined  in  detail  by 
comparing  model-generated  and  observed  river  runoff  and  precipitation  at  several 
locations. 

The  last  objective  of  this  thesis  is  to  examine  the  effects  of  doubling  the  current 
atmospheric  C02  level  and  to  determine  how  this  changes  monthly  river  runoff  for  wet, 
dry,  and  moderately  wet  river  basins  throughout  the  world,  particularly  the  impact  on 
snow  fall  and  snow  melt  for  high  latitude  rivers.  The  annual  results  of  Miller  and  Russel! 
[1992]  will  be  extended  to  examine  potential  seasonal  changes  in  future  runoff  due  to  the 
doubled  C02  climate  The  emphasis  will  be  on  the  high  latitude  river  basins  and  how 
spring  runoff,  due  to  snow  melt,  changes  in  the  doubled  C02  climate  and  how  the  different 
modeling  schemes  simulate  the  peak  runoff  season 

For  the  30  river  basins  in  this  study,  precipitation  and  river  runoff,  from  several 
different  simulations  will  be  examined.  The  emphasis  is  on  the  new  routing  scheme  and 
climate  change,  particularly  for  high  latitude  rivers.  The  results  are  analyzed  in  the  context 
of  several  other  studies  [ Gleick ,  1987;  Russell  and  Miller,  1990,  and  Kuhl  and  Miller, 
1992],  which  use  both  the  on-line  and  off-line  methods  to  calculate  river  runoff 
Suggested  improvements  on  how  the  model  generates  precipitation  and  river  runoff  are 
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Chapter  2.  Monthly  river  runoff  calculated  from  an  atmospheric  model 

2a.  Model  description  and  hydrology 

The  atmospheric  model  used  to  calculate  the  river  runoff  and  precipitation  was 
developed  at  the  NASA/Goddard  Institute  for  Space  Studies  and  is  described  by  Hansen 
et.  al,  [1983],  This  model  has  a  horizontal  resolution  of  4°  X  5°  in  latitude  and  longitude 
respectively,  and  a  vertical  resolution  of  nine  atmospheric  layers  The  model 
simultaneously  solves  the  equations  for  the  conservation  of  mass,  energy  and  momentum, 
and  the  equation  of  state  on  a  spherical  grid  Atmospheric  parameters  such  as  radiation  , 
atmospheric  gases  and  aerosols,  cloud  cover,  convection,  and  heat  and  momentum  are 
computed.  Other  parameters  include  ground  temperatures,  snow  depth  and  albedo,  and 
ground  hydrology.  The  primitive  equations,  conservation  of  momentum  (Newton's  2nd 
law  of  motion)  [Eqn  2.1],  conservation  of  mass  (continuity  equation)  [Eqn  2.2], 
conservation  of  energy  (1st  law  of  thermodynamics)  [Eqn  2.3],  and  the  ideal  gas  law  [Eqn 
2.4],  describing  the  state  and  motion  of  the  atmosphere  are  solved  numerically. 


—  =  -2Clx  V-p~'Vp  +  g  +  F 

(2.1) 

^  =  -pV.V  +  C-D 
dt 

(2.2) 

dl  dp 1  _ 

—  =  -p-t—+Q 

dt  Hdt 

(23) 

p~pRT 

(24) 
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where  V  is  the  velocity  relative  to  the  rotating  earth,  t  is  time,  Q  is  the  planet's  angular 
rotation  vector,  p  is  the  atmospheric  density,  g  is  the  apparent  gravity  [true  gravity  -  Q  x  ( 
Q  x  r)],  r  is  the  position  relative  to  the  planet's  center,  F  is  the  force  per  unit  mass,  C  is 
the  rate  of  creation  of  (gaseous)  atmosphere,  D  is  the  rate  of  destruction  of  atmosphere,  / 
is  the  internal  energy  per  unit  mass  [cvT\,  Q  is  the  heating  rate  per  unit  mass,  R  is  the  gas 
constant,  T  is  the  temperature,  and  cvFis  the  specific  heat  at  constant  volume. 

The  top  of  the  dynamical  atmosphere  is  fixed  at  lOmb.  In  the  model,  the  nine 
vertical  layers  include  two  in  the  boundary  layer,  five  in  the  troposphere,  and  two  in  the 
stratosphere.  At  the  surface,  grid  boxes  are  divided  into  land  and  ocean  fractions,  except 
for  the  model-simulation  C003,  in  which  grid  boxes  are  all  water  or  all  land  (see  Table 
3.1).  The  land  distribution  and  continental  topography  are  from  a  corrected  version  of 
Gates  and  Nelson  [1975],  Interactions  between  the  surface  and  atmosphere  are  computed 
separately  for  each  surface  type. 

The  atmospheric  model  separates  the  soil-moisture  storage  into  two  layers 
described  by  Hansen  et  al  [1983].  The  upper  layer  responds  immediately  to  evaporation 
and  precipitation  while  the  lower  layer  acts  as  a  seasonal  reservoir  The  rate  of  change  of 
moisture  in  the  upper  layer  is  given  by 

£KS  P-E-R  f  fV2-fV, 

dt  fx  t 

where  W,  (/  =  1  or  2),  is  the  ratio  of  available  water  to  the  field  capacity  (available  water  at 
saturation)  of  each  layer,  P  is  precipitation,  E  is  evaporation,  and  R  is  runoff  T  is  the  time 
constant  for  the  diffusion  of  moisture  between  the  two  soil  layers  and  is  the  field 
capacity.  The  rate  of  change  of  soil  moisture  storage  in  the  lower  layer  is  given  by 
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PW  f  W -W 

£^  =  A^l_£l  (26) 

dt  f2  r 

There  is  a  2-day  time  constant  for  diffusion  of  water  between  the  two  layers,  except 
during  the  growing  season  in  which  upward  diffusion  occurs  instantly  over  vegetated 
areas.  The  water  field  capacities  of  the  two  layers  depend  on  vegetation  characteristics 
taken  from  Matthews  [1983],  A  more  physically  realistic  formulation  of  the  surface 
hydrology  has  been  developed  by  Abramopoulos  et  al  [1988]  and  will  be  used  in  future 
simulations  with  the  GISS  model. 

The  computed  runoff  in  each  grid  box  depends  on  the  precipitation, 
evapotranspiration,  and  water  storage  within  the  land  portion  of  each  grid  box  The 
evapotranspiration  is  calculated  as  the  product 

E=ppCV(qs-qA)  (27) 

where  (3  is  a  dimensionless  efficiency  factor  for  evapotranspiration,  p  is  the  surface  air 
density,  C  is  a  dimensionless  drag  coefficient  that  depends  on  stability,  V  is  the  surface 
wind  speed,  qs  is  the  surface  saturation  specific  humidity  that  depends  on  ground 
temperature  and  the  surface  pressure,  and  qA  is  the  surface  air  specific  humidity  at  10m 
above  the  surface.  The  partition  of  rainwater  or  meltwater  between  the  fraction  that 
enters  and  remains  in  the  soil  and  the  fraction  moving  out  of  the  grid  box  depends  on  the 
soil  water  holding  capacity  and  the  soil  moisture.  The  water  leaving  the  grid  box  is  the 
runoff  [Miller,  1977],  Runoff./?  is  calculated  as 


R  =  maximum 


P  +  W- 


(2.8) 
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where  P  is  the  precipitation,  W  is  the  water  and  ice  in  the  first  layer,  and  Wc  is  the  water 
field  capacity.  The  factor  /?  in  (2.7)  is  assumed  to  equal  W  Wc  unless  the  ground  is  snow 
covered,  in  which  case,  P  =  1 .  The  coefficient  of  one-half  in  (2  8)  was  chosen  so  that  the 
computed  mean  annual  global  runoff  is  consistent  with  that  observed  [Hansen  et  al , 
1983],  In  the  parameterization  given  in  (2  8),  no  distinction  is  made  between  surface 
runoff  and  ground  water  runoff  that  leaves  each  grid  box.  Since  (2  8)  represents  the  total 
water  removed  from  a  grid  box,  it  combines  the  water  lost  from  the  two  components  of 
runoff  but  without  any  time  lag,  thus,  runoff  occurs  during  precipitation  and  stops  when 
the  precipitation  stops.  Runoff  is  a  continuous  function  of  precipitation.  In  the  model 
precipitation  can  be  in  the  form  of  rain  or  snow.  The  snow  depth  for  any  given  area  is 
computed  as  the  balance  of  snowfall,  melting  and  sublimation  [Hansen  et  al,  1983]  If  the 
air  temperature  is  less  than  0°C,  precipitation  falls  as  snow.  If  the  temperature  of  the 
upper  layer  of  soil  is  less  than  or  equal  to  0°C,  the  srn  .  depth  will  increase,  otherwise  the 
snow  melts  and  the  upper  soil  temperature  decreases  The  snow  melt  is  then  included  in 
the  runoff. 

The  drainage  basins  for  this  study  were  defined  by  Russell  and  Miller  [1990]  using 
the  maps  of  Korzoun  et  al  [1977]  and  the  Times  Atlas  of  the  World  [1967]  The  model's 
4°  X  5°  grid  boxes  were  divided  into  four  2.5°  X  2°  grid  boxes  with  the  total  area  of  the 
river  drainage  basin  equal  to  the  sum  of  the  areas  of  the  2.5°  X  2°  grid  boxes.  If  a 
particular  grid  box  was  assigned  to  a  river  basin,  all  the  runoff  from  that  grid  box  was 
assigned  to  the  river  flow  In  Kuhl  and  Miller  [1992],  computed  monthly  runoff  was 
defined  as  the  sum  of  the  runoff  from  each  grid  box  within  the  river  basin  during  the 
month.  This  is  not  the  same  as  observed  runoff,  which  is  defined  as  the  monthly  water 
flow  at  the  mouth  of  the  river.  Comparisons  between  model-generated  runoff  and 
observed  runoff  are  discussed  in  Miller  and  Russell  [  1 992]  and  Kuhl  and  Miller  [  1 992] . 
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2b.  River  routing  scheme. 

A  river  routing  scheme  has  been  developed  to  simulate  water  flow  from  one  grid 
box  to  another.  This  allows  a  time  lag  in  the  routing  of  river  runoff,  in  particular  it  allows 
one  to  calculate  river  runoff  at  the  mouth  of  the  river  basin  and  more  easily  compare  it 
with  the  observed  data  The  calculation  of  the  monthly  runoff  in  the  routing  scheme 
requires  directed  paths  connecting  grid  boxes  within  the  river  basin  to  grid  boxes  which 
are  successfully  closer  to  the  river's  mouth.  This  is  discussed  in  Miller  et  al.  [1992],  The 
rate  at  which  water  moves  between  grid  boxes  depends  on  many  factors  such  as  the  slope, 
the  volume  of  the  water  in  the  river,  bottom  composition,  shape,  width,  depth,  and 
distanced  traveled,  of  which  only  slope,  river  volume,  and  distanced  traveled  are  included 
in  the  model  simulations. 

River  direction  files  were  created  using  the  2  5°  X  2°  resolution  described  in  the 
previous  section.  Each  grid  box  contains  a  value  from  0  through  8  defining  the  direction 
of  flow  within  a  grid  box.  A  value  of  zero  means  no  flow  out  of  the  grid  box  Values 
greater  than  zero  means  water  flows  out  of  the  grid  box  via  other  rivers  to  an  adjacent 
grid  box  downstream.  A  value  of  one  represents  water  flow  downstream  to  the  northeast, 
and  values  from  2  through  8  represent  water  flow  downstream  to  the  north,  northwest, 
west,  southwest,  south,  southeast,  and  east,  respectively.  Direction  files  were  extracted 
from  world  maps  [Korzoun  et  al,  1977], 

Once  the  direction  files  have  been  established,  the  next  step  is  to  determine  the 
volume  of  the  river  water  in  each  grid  box  and  the  flow  rate  between  each  grid  box.  This 
approach  is  modeled  after  continuous  streamflow  simulation  models  discussed  by  Singh 
[1989],  The  volume  of  water  in  a  grid  box  after  a  time  step  At  is  given  by 


v(t  +  A ,)  =  F(l)  +(/-,  -  Fj A I  +  K,R,&  <  +  '  <2  9> 
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where  Fin  and  Fout  are  the  flow  rates  into  and  out  of  the  grid  box,  Ks  and  Kg  are  rate 
constants  for  surface  and  ground  water  runoff,  respectively,  and  Rs  and  Rg  are  reservoirs 
of  surface  and  groundwater  runoff  For  each  grid  box,  flow  into  the  box  can  be  from  all 
directions,  however,  flow  out  of  a  grid  box  can  only  be  in  one  direction  and  is  assumed  to 
be  of  the  form 

F,  =£!,/(/>,  Ay/ As,  (2  10) 

where  /;  is  the  slope  of  the  grid  box,/is  a  function  of  ip  Q,  (s*1)  is  a  rate  coefficient,  Ay  is 
the  north-south  width  of  the  grid  box,  and  Asy  is  the  distanced  traveled  by  the  river  across 
the  grid  box.  The  subscripts  indicate  that  the  variables  may  differ  from  grid  box  to  grid 
box 

The  rate  coefficient,  fij„  depends  on  soil  type  and  depth,  characteristics  of  the 
drainage  basin  network,  and  vegetation  cover,  however,  will  be  taken  as  a  constant 
independent  of  these  variables.  Equation  2 .10  is  based  on  the  Muskingum  method 
[Linsley  et  al ,  1982]  and  represents  a  simplified  form  of  the  runoff  routing  scheme  used  by 
Vorosmarty  et  al  [1989],  As  discussed  in  Miller  et  al  [1992],  Q,  is  dependent  on  grid 
resolution.  If  the  grid  size  increases  than  Qj  will  decrease.  Miller  et  al  [1992]  have  given 
an  alternative  formulation  to  Eqn.  2.10  in  which  they  define  a  basinwide  turnover  rate 
which  is  independent  of  grid  resolution. 

The  slope,  ip  is  dimensionless  and  is  calculated  as  the  ratio  of  the  height  difference 
between  the  upstream  and  downstream  grid  boxes  and  the  distanced  traveled.  The 
distanced  traveled  depends  on  the  direction  of  flow  in  a  grid  box,  east-west,  north-south, 
or  diagonally,  which  is  found  on  a  directional  map  file  discussed  before.  If  the  river  flows 
in  the  north-south  direction,  then  Asy  =  Ay,  and  equation  2  10  reduces  to  F}  =  Vp  If 
the  flow  is  east-west,  than  the  distanced  traveled  will  decrease  as  one  moves  northward 
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due  to  the  convening  meridians.  If  the  river  flows  diagonally,  then  i\s2  =  {Ax2  +  Ay2) 
Miller  and  Russell  [1992]  discuss  the  effects  of  grid  box  size  and  the  direction  and 
distance  traveled. 

The  next  step  in  the  routing  scheme  is  to  include  topography.  Grid  boxes  with 
steeper  channel  slopes  are  assumed  to  flow  faster.  Variations  in  topography  are  shown  in 
Miller  et  al  [1992]  for  four  of  the  rivers  discussed  in  this  paper.  In  this  work 

/(>',)  =  (2.11) 


where  imin  is  the  minimum  slope  allowed  for  a  grid  box,  hence,  imin  <  /;  <  imax  and  imin  = 
20m/ As  and  imax  =  200m/ As, 

The  routing  scheme  described  in  equations  2  9  and  2.10  allows  river  water  to 
move  more  realistically  from  its  origin  to  the  mouth  of  the  river  basin  It  also  allows  one 
to  calculate  the  runoff  at  any  location  on  the  river,  the  total  volume  of  water  in  each  grid 
box  and  in  the  entire  river  basin,  and  the  total  instantaneous  runoff  into  the  river  Kuhl 
and  Miller  [1992]  examined  monthly  runoff  variation  of  the  total  runoff  (sum  over  all  grid 
boxes)  for  16  river  basins.  The  next  chapter  will  extend  their  work  and  compare  runoff  at 
the  river  mouth  with  the  total  runoff. 
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Chapter  3.  Monthly  River  RunoiT  and  Precipitation 


In  this  chapter,  monthly  river  runoff  and  precipitation  generated  by  several 
different  simulations  with  the  GISS  GCM  are  compared  for  several  major  liver  basins 
Table  3. 1  shows  how  the  horizontal  grid,  sea  surface  temperature,  and  run  length  for  the 
different  simulations  The  purpose  of  the  comparison  is  to  provide  some  insight  into  the 
sensitivity  and  variability  of  the  model 


Table  3.1.  Comparison  of  the  different  model  simulations  used  in  this  paper.  For  simplicity,  model 
simulations  will  be  referred  to  bv  the  model  number  show  n  in  the  first  column. 


Model 

Horizontal 

Grid1 

CO 2  Level  j 

Simulation 

Length 

Climatological 
Sea  surface  temps  | 

Fractional  Grid 

Present 

Doubled 

Yes 

No 

Yes 

No 

848 

B 

■1 

■/ 

✓ 

B100 

B 

✓ 

mm 

✓ 

C003 

C 

✓ 

_  5yr 

✓ 

✓2 

A51 

B 

mm 

✓ 

947 

B 

V 

mm 

✓3 

Horizontal  grid  system  used  based  on  Hansen  et  al.  [1983] 

2  Non-fractional  grid  Each  gnd  box  is  either  all  water  or  all  land 

3  Sea  surface  temperatures  are  interpolated  from  equilibrium  simulations  at  8°  X  10°  resolution  in  which  SST  were 
predicted  in  Hansen  et  al.  [  1 984], 


Several  different  three-year  and  five-year  simulations  are  available  with  somewhat 
different  model  formulations  for  the  present  climate  Because  of  the  relatively  short 
simulations,  it's  not  certain  that  the  models  represent  long  term  climatology  for  runoff  and 
precipitation.  If  all  the  simulations  are  in  agreement  with  each  other,  than  runoff  and 
precipitation  for  the  present  climate  should  be  representative  of  the  long  term  climatology 
However,  if  the  simulations  differ,  there  would  be  less  confidence  in  the  results  for  that 
basin. 
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3a.  Observed  river  runoff  and  precipitation 

Figure  3.1  shows  the  location  of  the  mouth  of  all  rivers  used  for  this  study  The 
observed  monthly  river  runoff  in  this  study  is  the  flow  at  the  location  nearest  the  mouth 
within  the  river  basin  and  is  based  on  UNESCO  data.  During  the  1960’s,  the  thirteenth 
General  Conference  of  UNESCO  launched  the  International  Hydrological  Decade,  (IHD) 
1965-1974.  Part  of  the  program  focused  on  recording  and  collecting  rates  of  river  runoff 
at  various  locations  within  the  drainage  basin.  The  number  of  years  used  for  river  runoff 
averages  vary  between  river  basins.  The  majority  of  the  river  basins  have  more  than  20 
years  of  data,  while  the  remaining  basins  contain  only  four  years  of  data  The  observed 
precipitation  accumulated  over  the  model's  drainage  basin  area  is  taken  from  two  studies. 
Shea  [1986]  and  Legates  and  Willmott  [1990], 

In  previous  studies  of  Rime//  and  Miller  [1990],  Kuhl  [1990],  Miller  and  Russell 
[1992],  and  Kuhl  and  Miller  [1992],  the  observed  precipitation  was  based  on  Shea 
[1986],  Figures  3.2-3  5  show  comparisons  of  observed  monthly  precipitation  between 
Legates  and  Willmott  [1990]  and  Shea  [1986]  for  several  of  the  world's  major  river  basins 
For  the  majority  of  the  river  basins,  annual  precipitation  between  the  two  studies  is  within 
10%  with  Legates  and  Willmott  [1990]  averaging  slightly  higher  annua!  amounts  The 
greatest  disparity  between  the  two  observed  studies  occurs  in  the  wet  climates  (Figure 
3.2),  however,  despite  an  increase  in  precipitation  amounts  of  30%  for  several  of  the 
rivers,  the  overall  monthly  variation  in  precipitation  is  similar  for  Legates  and  Willmott 
[1990]  and  Shea  [1986],  This  is  true  for  all  the  rivers  in  Figures  3.2-3  5,  thus,  observed 
precipitation  in  the  remainder  of  this  study  was  taken  from  Legates  and  Willmott  [1990], 
which  is  the  more  recent  of  the  two  studies  and  contains  more  observed  data 
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3b.  River  runoff  for  the  world's  major  rivers 

Modei-generated  monthly  river  runoff,  R,ol,  discussed  in  Kuhl  and  Miller  [1992] 
was  calculated  by  summing  the  runoff  from  each  grid  box  over  the  river  basin,  as 
discussed  in  the  previous  chapter.  This  is  different  from  observed  runoff  at  the  mouth  of 
the  river  basin,  R^,  which  represents  the  flow  of  water  at  the  mouth  after  the  water  has 
moved  through  the  entire  river  basin,  thus  R,o1  should  lead  A  time  dependent  river 
routing  scheme  has  been  developed  by  Miller  el  at.  [1992]  to  allow  the  calculation  of 
model-generated  river  runoff  anywhere  in  the  river  basin,  as  discussed  in  chapter  2  To 
better  understand  how  the  new  routing  scheme  works,  the  B100  model  was  used  to 
calculate  R^  and  R,ot  Figures  3  6-3  9  show  model-generated  and  observed  precipitation 
and  runoff  for  several  of  the  world’s  major  river  basins  The  remainder  of  this  paper  will 
refer  to  river  runoff  at  the  mouth  as  R,,,  and  the  sum  of  river  runoff  basinwide  as  R,0I 

The  timing  of  the  model-generated  runoff  without  the  routing  scheme,  Rlot.  is 
primarily  based  on  the  model's  precipitation  As  explained  in  chapter  2,  runoff  continues 
as  long  as  the  precipitation  continues  and  when  precipitation  stops  runoff  stops.  A  good 
example  of  how  the  model  acts  without  the  routing  scheme  is  in  the  wet  climates  (Figure 
3  6)  where  there  is  no  snowfall  Figure  3  6a  shows  that  the  variation  of  the  precipitation 
for  the  Congo  River  is  similar  to  that  of  Rtol,  although  precipitation  is  too  large  Figure 
3  6a  shows  that  the  maximum  model-generated  precipitation  occurs  in  March  and 
minimum  model-generated  precipitation  occurs  in  July  The  model-generated  runoff,  R,ot. 
follows  the  same  pattern.  However,  for  the  model-generated  runoff,  R,,,,  the  river  runoff 
in  March  is  lower  and  the  runoff  from  April  through  July  is  higher  at  the  mouth  which 
indicates  that  the  water  contributing  to  the  March  peak  in  R,o1  reaches  the  mouth  in  the 
following  months  The  opposite  is  true  for  the  period  August  through  October,  the  period 
between  the  minimum  and  next  maximum  precipitation,  in  which  R,pt  follows  the  same 
pattern  as  the  precipitation  curve,  however,  R^  shows  the  river  runoff  lagging  behind 
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slowly  building  up  for  a  maximum  in  January  The  model-generated  river  runoff,  Rm, 
follows  the  general  patterns  of  the  observed  runoff  Similar  results  are  shown  in  Figuies 
3.6b  through  3.8  for  the  Mekong  and  other  rivers 

The  effect  of  the  new  routing  scheme  is  illustrated  best  in  the  higher  latitude  river 
basins  where  snow  occurs  and  snow  melt  is  a  major  contributor  to  river  runoff.  Where 
snow  accumulates  during  the  winter,  river  runoff  maxima  usually  occur  when  the  snow 
melts  during  the  spring  and  early  summer.  Figure  3.9  shows  monthly  precipitation  and 
runoff  from  the  Yenesei  and  Amur  river  basins  located  in  northeastern  Siberia.  Although 
the  annual  variation  of  the  model-generated  precipitation  is  in  good  agreement  with  the 
observed,  for  both  rivers,  it  is  too  high  in  the  spring  Precipitation  peaks  in  July  the 
Yenesei  Rivers,  however,  unlike  the  wet  climates  discussed  before,  the  runoff  maximum 
occurs  prior  to  the  precipitation  maximum.  This  is  due  to  the  snow  melt  in  spring  and 
early  summer  In  the  model  the  maximum  snow  melt  occurs  in  March  and  April,  with 
model-generated  runoff  peaks  occurring  during  the  same  months  for  R,ot  Results  are 
similar  for  the  other  high  latitude  rivers.  With  the  routing  scheme  in  place,  model¬ 
generated  runoff  is  allowed  to  flow  downstream  towards  the  mouth  of  the  river 
throughout  the  snow  melt  season  peaking  in  late  spring  and  1^,  is  in  reasonably  good 
agreement  with  the  observed 

The  new  routing  scheme  simulates  water  moving  downstream  and  enables 
modeler's  to  simulate  the  magnitude  and  timing  of  the  runoff  peak  The  new  routing 
scheme  also  allows  modeler's  to  choose  any  location  within  the  river  basin  and  analyze  the 
instantaneous  runoff  flow  at  that  specific  location.  The  routing  scheme  does  not  deal 
directly  with  movement  of  water  within  each  grid  box.  This  is  done  using  the  grid  box 
hydrology  scheme  discussed  in  chapter  2.  It  should  be  noted  that  although  the  new 
routing  scheme  produces  changes  in  the  timing  of  the  runoff  peak,  only  monthly  data  exist 
If  daily  data  were  available,  the  timing  of  runoff  peaks  would  be  further  resolved  giving 
modeler's  a  more  precise  tool  for  forecasting  the  magnitude  and  timing  of  river  runoff 
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Also  of  note  is  that  the  routing  coefficients  are  constant  and  have  not  been  tuned  for 
individual  river  basins.  There  is  some  variation  of  flow  rates  based  on  topography 
Optimal  results  may  not  be  possible  until  the  routing  coefficients  are  tuned  to  individual 
basins. 

The  different  model  simulations  for  the  present  climate  used  in  this  study  were 
listed  in  Table  3.1.  Table  3.2  shows  model-generated  precipitation  and  runoff  for  30  river 
basins  comparing  the  different  model  simulations  used  in  this  paper  along  with  the 
observed  precipitation  and  runoff.  As  discussed  before,  longer  simulations  would  be 
better  for  comparing  model-generated  river  runoff  with  observed  river  runoff.  Overall,  all 
four  model  simulations  produce  too  much  precipitation  and  are  in  general  agreement  with 
each  other  more  than  with  the  observed  precipitation  With  the  exception  of  the  wet 
climate  river  basins,  river  runoff  is  too  high  for  all  four  simulations 

The  four  simulations  use  two  of  the  three  horizontal  grid  schemes  described  in 
Hansen  et  a!.,  [1983],  scheme  B  and  scheme  C.  In  the  B  and  C  grid  schemes,  the  pressure 
gradient  and  velocity  divergence  are  computed  over  AX,  resulting  in  more  accurate 
representation  of  the  geostrophic  adjustment  [Arakawa,  1972],  The  difference  between  he 
B  and  C  scheme  is  that  the  winds  are  directly  computed  on  grid  C,  but  on  grid  B  it  is 
necessary  to  average  the  winds.  Less  averaging  of  the  winds  should  lead  to  more  accurate 
transports  of  heat  and  water  with  the  C  grid  scheme.  More  detailed  explanations  on  the 
different  grid  schemes  are  in  Hansen  et  al.  [1983], 

Figures  3.10  through  3.13  show  comparisons  of  model -generated  precipitation  and 
runoff  for  four  river  basins.  Model  generated  runoff  shown  is  R,ol  One  interesting  case  is 
the  Indus  River  (Fig  3.12)  in  which  the  model  C003  is  the  only  one  to  show  a  summer 
maximum  in  precipitation  due  to  the  monsoon  season.  The  three  other  models  have 
precipitation  maxima  in  the  winter  months.  Model  C003  is  the  only  model  simulation  to 
use  a  non-fractional  grid  Each  grid  box  is  either  water  or  land  Model  C003  also  uses  the 
C  horizontal  grid  scheme  discussed  above.  These  two  factors  may  cause  the  difference  in 
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the  monthly  precipitation  in  the  Indus  River  basin  Overall,  the  annual  model-generated 
precipitation  and  river  runoff  for  all  four  simulations  are  within  8%  and  1 3%,  respectively, 
of  each  other.  The  monthly  precipitation  and  river  runoff  are  for  the  most  part  in 
agreement  with  one  another.  The  only  difference  is  the  amount  of  precipitation  and  river 
runoff  flow.  The  largest  discrepancies  occur  in  the  dry  basins,  as  seen  with  the  Indus 
River,  in  which  the  four  simulations  differ  in  precipitation  and  river  runoff  by  50%  and 
235%,  respectively. 

Throughout  this  paper  comparisons  of  annual  and  monthly  precipitation  and  runoff 
are  discussed.  Annual  model-generated  precipitation  and  runoff  may  be  within  5%  of  the 
annual  observed  precipitation  and  runoff,  as  in  the  case  of  the  Mississippi  River  (within 
2%  of  the  observed  precipitation)  and  Lena  River  (within  4%  of  the  observed  runoff), 
however,  monthly  precipitation  and  river  runoff  could  be  incorrect  as  shown  for  the 
Mississippi  and  Lena  rivers  (Fig  3.14)  The  Lena  River,  though  within  4%  of  the  annual 
observed  river  runoff,  peaks  to  soon  in  the  model-generated  river  runoff,  R^  The 
Mississippi  River  has  too  much  model-generated  precipitation  in  the  spring  and  too  little 
model-generated  precipitation  in  the  fall.  These  are  examples  of  cases  in  which  annual 
precipitation  and  runoff  are  simulated  well,  but  monthly  precipitation  and  runoff  are  not  as 
well  simulated. 

The  factor  Cl  in  Eqn.  2.10  is  a  key  variable  in  calculating  the  runoff  at  a  rivers 
mouth.  In  that  equation  Q  w^s  taken  to  be  constant  globally  As  Q  continues  to  increase, 
and  if  the  time  step  is  relatively  short,  the  variation  of  runoff  at  the  mouth  will  approach 
that  of  Rtot.  As  O  gets  small,  water  will  build  up  in  the  grid  box  smoothing  the  flow  rate 
over  time.  The  runoff  rates  given  by  Eqn.  2.10  for  individual  grid  boxes  could  be  modified 
further  by  allowing  Cl  to  vary  depending  on  the  physical  characteristics  of  a  particular  grid 
box.  These  might  include  soil  type  and  depth,  vegetation  type,  and  topography  or  other 
characteristics  of  the  river  basin 
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Another  problem  which  may  exist  in  analyzing  individual  river  basins  is  that  the 
precipitation  and  runoff  are  occurring  in  the  correct  part  of  the  river  basin  With  the  new 
routing  scheme,  it  is  possible  to  examine  runoff  for  separate  components  of  a  larger  river 
basin.  For  example,  the  Mississippi  River  can  be  divided  into  its  components,  the 
Missouri,  Ohio,  Illinois,  and  Arkansas  Rivers  for  further  analysis.  The  new  routing 
scheme  allows  us  to  choose  any  location  within  the  basin  and  simulate  runoff  there  The 
next  chapter  will  discuss  the  precipitation  and  runoff  for  the  Mississippi  River  and  its 
components. 


Table  3.2.  Annual  observed  and  model-generated  runofT  (Km^/yr)  and  precipitation  (Km3/yr)  for 
the  world's  major  rivers  for  several  different  model  simulations. 


Climate  Type 
River 

OBS 

Precipitation 

848  B100  C003 

AS1 

OBS 

848 

Runoff 

B100 

COOL 

ASJ 

Wet 

Amazon 

13942 

12296 

12842 

13241 

14224 

4886 

2338 

2479 

4055 

3079 

Congo 

5516 

8841 

7635 

6446 

8470 

1398 

2170 

1574 

1643 

2094 

Orinoco 

2193 

2173 

2221 

2040 

2005 

794 

475 

615 

707 

466 

Mekong 

1477 

1996 

2112 

1758 

2133 

449 

712 

/U9 

410 

710 

Magdalena 

758 

863 

1176 

1095 

1056 

213 

314 

520 

581 

304 

Sao  Francisco 

799 

1160 

987 

821 

1273 

83 

211 

178 

189 

272 

Average 

4114 

4555 

4496 

4234 

4861 

1304 

1037 

1013 

1264 

1154 

Moderately  Wet 

Yangtze 

1985 

3436 

3167 

2787 

3465 

792 

1304 

1239 

898 

1499 

Mississippi 

3039 

2645 

2965 

2127 

2981 

498 

517 

646 

395 

660 

LaPlata  (Parana) 

3662 

2955 

2629 

3905 

3164 

— 

(407) 

(369) 

(1444) 

(474) 

St  Lawrence 

1070 

1204 

1144 

1053 

1102 

214 

462 

361 

505 

322 

Danube 

1192 

995 

1349 

1289 

1683 

— 

(298) 

(352) 

(309) 

(516) 

Columbia 

507 

753 

660 

669 

774 

172 

304 

227 

258 

248 

Zambesi 

1218 

1600 

1319 

1058 

1356 

105 

256 

177 

204 

198 

Fraser 

179 

288 

278 

362 

307 

87 

159 

127 

215 

134 

Nile 

1960 

3508 

3601 

3228 

3659 

— 

(587) 

(634) 

(758) 

(656) 

Niger 

1618 

1917 

2306 

2136 

2116 

(350) 

(403) 

(418) 

(331) 

Average 

1643 

1930 

1942 

1861 

2061 

311 

500 

463 

413 

510 

Dry 

Indus 

392 

660 

314 

543 

592 

76 

302 

116 

198 

204 

Tigns-Euphrates 

497 

500 

503 

442 

497 

— 

179) 

(86) 

(126) 

(75) 

Yellow 

541 

1403 

1347 

1088 

1454 

— 

(514) 

(497) 

(275) 

(563) 

Colorado 

205 

418 

408 

198 

339 

12 

83 

87 

33 

70 

Murray 

528 

727 

639 

336 

585 

8 

117 

98 

44 

97 

Average 

433 

742 

642 

521 

693 

32 

219 

177 

135 

202 

Hieh  Latitude  Dry 

Yenesei 

1126 

1480 

1646 

1788 

1371 

558 

499 

586 

710 

461 

Lena 

926 

1407 

1590 

1839 

1486 

516 

540 

638 

779 

535 

Ob 

1354 

1248 

1131 

1440 

1066 

388 

503 

455 

605 

464 

Amur 

1120 

1364 

1484 

1635 

1394 

309 

316 

357 

475 

273 

Mackenzie 

607 

1172 

1159 

1134 

1242 

264 

560 

514 

599 

644 

Yukon 

333 

779 

707 

886 

715 

197 

492 

519 

713 

510 

Sevemay  Dvina 

206 

224 

237 

284 

241 

107 

118 

117 

152 

132 

Kolyma 

258 

470 

559 

514 

506 

71 

337 

434 

411 

366 

Indigirka 

100 

205 

233 

214 

228 

49 

119 

146 

135 

134 

Average 

670 

928 

972 

1082 

917 

273 

387 

418 

509 

391 

Total  Average 

1644 

1956 

1945 

1879 

2050 

510 

550 

538 

621 

578 

Observed  precipitation  is  from  Legates  and  Willmott  (1990).  Observed  runofT  is  from  UNESCO  (1969,1974,1985) 
{ — )  Indicates  river  basins  in  which  observed  runoff  station  is  not  near  the  mouth  of  the  river.  Values  in  (  )  are 
not  included  in  climate  type  and  total  runofT  averages 
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AMAZON _  CONGO  ORINOCO 


21 


Figure  3  2  Observed  monthly  precipitation  from  Legates  and  Willmott  [1990]  and  Shea  [1986]  for  the  (a)  Amazon 
Orinoco,  (d)  Mekong,  (e)  Magdalena,  and  (f)  Sao  Francisco  rivers. 


YANGTZE  ZAMBESI  MISSISSIPPI 


Mississippi,  (d)  Nile,  (e)  LaPlata  (Parana),  and  (0  Niger  rivers. 
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MACKENZIE  SEVERNAY 
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Figure  3.5  Observed  monthly  precipitation  from  Legates  and  Willmott  [1990]  and  Shea  [1986]  for  the  (a)  Mackenzie,  (b)  Sevemay 
Dvina,  (c)  Ob,  (d)  Amur,  (e)  Lena,  and  (0  Yenesei  rivers. 


Precipitation  Runoff 


Figure  3.8.  Model-generated  (M)  and  observed  (O)  precipitation  and  river  runoff  for  the  (a)  Indus  and  (b)  Yellow  rivers  R,0(  is  model 
river  runoff  averaged  over  the  entire  river  basin  while  Rm  is  model  river  runoff  at  the  mouth  using  the  new  routing  scheme.  Observed 
precipitation  is  from  Legates  and  Willmott  [1990]  and  observed  runoff  is  from  UNESCO  [1969,  1974,  1985], 
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Figure  3.11.  Model-generated  (A51,  B100,  C003,  848)  and  observed  (O)  precipitation  and  river  runoff  for  the  Mississippi  River. 
Observed  precipitation  is  from  Legates  and  Willmott  [1990]  and  observed  runoff  is  from  UNESCO  [1969,  1974,  1985],  A51,  B100, 
C003,  and  848  refer  to  the  simulations  shown  in  Table  3.1.  River  runoff,  Rtot,  is  averaged  over  the  entire  river  basin 


INDUS  INDUS 
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Figure  3.12.  Model-generated  (A51,  B100,  C003,  848)  and  observed  (O)  precipitation  and  river  runoff  for  the  Indus  River  Observed 
precipitation  is  from  Legates  and  Willmott  [1990]  and  observed  runoff  is  from  UNESCO  [1969,  1974,  1985],  A51,  B100,  C003,  and 
848  refer  to  the  simulations  shown  in  Table  3. 1 .  River  runoff,  Rt0(,  is  averaged  over  the  entire  river  basin. 
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Chapter  4.  Mississippi  River  Basin 

Hydrologists  have  generally  focused  on  calculating  river  runoff  at  smaller  scales 
(about  1  kilometer)  than  those  considered  by  global  atmospheric  modelers  Global 
atmospheric  modelers  have  used  GCMs  to  calculate  river  runoff  at  much  larger  scales. 
Vorosmarty  et  al.  [1989]  have  developed  a  continental  drainage  basin  runoff  scheme  and 
applied  it  to  the  Amazon  River  basin  using  observed  data  as  input.  The  model  calculates 
water  flow  from  one  grid  box  to  another,  however,  it  is  based  on  observed  data  and  can 
not  be  used  directly  to  calculate  potential  changes  in  river  runoff  due  to  global  warming 
The  new  routing  scheme  discussed  in  chapters  2  and  3  allows  the  calculation  of  runoff  at 
any  grid  box  in  a  river  basin.  The  new  routing  scheme  was  used  in  the  analysis  of  the 
Mississippi  river  and  its  tributaries. 

Figure  4.1  shows  a  map  of  the  Mississippi  River  basin  divided  into  2°  X  2.5°  grid 
boxes  and  the  direction  of  the  downstream  flow  for  each  grid  box.  The  map  is  also 
divided  into  the  major  tributaries  of  the  Mississippi  River,  the  Missouri,  Arkansas,  Illinois, 
and  Ohio  rivers.  For  each  river  basin,  runoff  moves  from  grid  box  to  grid  box  within  the 
basin  and  is  calculated  at  the  grid  box  or  sum  of  grid  boxes  which  leave  the  basin  and  enter 
into  the  Mississippi  River.  Observed  and  model-generated  precipitation  and  snow  mass 
are  interpolated  from  a  4°  X  5°  grid  box  into  four  2°  X  2.5°  grid  boxes.  Precipitation  and 
snow  mass  were  assumed  to  be  equal  throughout  the  grid  box.  Hence,  if  the  precipitation 
averaged  2  mm  per  day  in  a  4°  X  5°  grid  box  than  it  would  average  the  same  2  mm  per 
day  in  each  2°  X  2.5°  grid  box. 

The  drainage  area  of  each  basin  is  the  sum  of  the  areas  of  all  grid  boxes  within  the 
river  basin.  The  area  of  each  grid  box  was  calculated  as  follows: 


Area  =(2  x2.5)  x(l  10£m)‘  x  cos# 


(4  1) 
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where  0  is  the  latitude  of  the  northern  edge  of  a  grid  box  Table  4  1  shows  model  and 
observed  areas  of  each  river  basin  and  three  locations  within  the  Mississippi  River  basin, 
Keokuk,  Iowa,  St  Louis,  Missouri,  and  Vicksburg,  Mississippi  The  model  to  observed 
ratio  is  near  one  for  all  except  the  three  smallest  basins,  Arkansas,  Illinois,  and  Upper 
Mississippi  rivers. 


Table  4.1.  Model  and  observed  drainage  basin  areas  (Km2)  for  the  Mississippi  River  and  it's 
tributaries.  The  last  column  is  the  ratio  of  the  model  to  the  observed  drainage  area.  _ 


River 

Area 

Area 

Location 

Model 

Observed 

Ratio 

Missouri  /  Boonville,  Missouri 

1283025 

1299403 

099 

Arkansas  /  Dardanelle,  Arkansas 

428955 

398005 

1  08 

Ohio  /  Metropolis,  Illinois 

518818 

525770 

099 

Illinois  /  Meredosia,  Illinois 

88480 

67412 

1.31 

Upper  Mississippi  /  Keokuk,  Iowa 

385552 

308210 

1.25 

Central  Mississippi  /  St  Louis,  Missouri 

1803400 

1805230 

1.00 

Lower  Mississippi  /  Vicksburg,  Mississippi 

2991969 

2953895 

1.01 

Observed  runoff  is  from  the  US  Geological  Survey  [  1 987] 

The  Illinois  and  Missouri  rivers  enter  the  Mississippi  River  before  reaching  St 
Louis  and  the  Ohio  and  Arkansas  rivers  enter  the  Mississippi  River  between  St  Louis  and 
Vicksburg.  Figures  4.2-4  4  show  observed  and  model-generated  precipitation  for  the 
various  river  basins  within  the  Mississippi.  The  figure  for  the  lower  Mississippi  river  at 
Vicksburg  (Fig  4.4)  is  the  same  as  Fig  3.7a.  As  discussed  before,  the  overall  annual 
model-generated  precipitation  is  within  5%  of  the  observed  precipitation,  however,  there 
is  too  much  precipitation  in  the  spring  and  too  little  in  the  summer  and  fall  The  reason  to 
further  divide  the  Mississippi  River  into  it's  tributaries  is  to  find  the  origin  of  the 
precipitation  and  runoff  and  how  it  compares  to  observed  values. 

In  general,  the  model  is  too  dry  in  the  eastern  part  and  too  wet  in  the  western  part 
of  the  Mississippi  basin.  The  high  precipitation  in  the  spring  is  caused  by  too  much 
precipitation  being  generated  in  the  Missouri  River  basin,  which  encompasses  44%  of  the 
Mississippi  River  basin.  The  late  summer  and  fall  differences  are  caused  by  too  little 
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precipitation  being  generated  in  the  upper  Mississippi,  Arkansas  and  Ohio  rivers,  42%  of 
the  Mississippi  River  basin.  Therefore,  although  the  annual  precipitation  is  within  5%  of 
the  observed,  the  monthly  precipitation  generated  by  the  mode!  is  occurring  in  the  wrong 
regions  of  the  basin.  This  also  affects  the  monthly  river  runoff  at  the  mouth  The 
observed  river  runoff  peaks  in  March  for  the  Ohio  River  and  April  for  the  Mississippi 
River  at  Vicksburg  while  the  model  shows  river  runoff  peaks  in  May  at  Vicksburg, 
following  an  April  peak  for  the  Missouri  River.  Although  the  observed  and  model- 
generated  river  runoff  peaks  for  Vicksburg  are  off  by  one  month,  the  model  is  able  to 
simulate  the  lag  as  water  moves  downstream  to  the  mouth  of  the  river  basin 

The  model  generates  too  much  annual  river  runoff  (16%)  for  the  entire  river  basin 
Figures  4.2a,  4  3a,  and  4.4  show  the  upper,  central,  and  lower  stations  along  the 
Mississippi  River.  The  model  generates  10%,  170%,  and  16%,  too  much  river  runoff, 
respectively,  for  all  three  stations.  The  model-generated  runoff  is  too  high  for  the 
Missouri  River  and  too  low  for  the  Ohio  River.  The  problem  with  the  Missouri  River  is 
that  the  model  generates  too  much  precipitation  from  October  through  July,  some  of 
which  is  in  the  form  of  snow.  The  model  generates  too  little  precipitation  from  July 
through  February  for  the  Ohio  River  causing  the  low  model -generated  runoff  The 
combination  of  too  much  runoff  for  the  Missouri  River  and  too  little  runoff  for  the  Ohio 
River  causes  the  annual  runoff  to  be  within  1 6%  of  the  observed  runoff  near  Vicksburg.  A 
separate  analysis  of  the  Mississippi's  tributaries  can  show  why  monthly  and  annual 
simulations  differ  in  a  particular  river  basin. 

A  useful  toot  in  examining  river  runoff  is  comparing  the  annual  runoff  coefficient, 
which  is  defined  as  the  ratio  of  annual  runoff  to  annual  precipitation,  for  observed  and 
model-generated  river  runoff.  The  observed  and  model  runoff  coefficients  for  the 
Missouri  River  are  .08  and  .29,  respectively.  The  observed  and  model  coefficients  for  the 
Ohio  River  are  .38  and  .17,  respectively.  The  significantly  different  observed  runoff 
coefficients  for  the  Missouri  and  Ohio  rivers,  which  is  not  so  prominent  in  the  model, 
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indicates  that  the  runoff  generation  for  the  two  basins  given  by  Eqn  2.8  should  be 
examined  further.  Possible  differences  between  the  two  basins  are  that  most  of  the 
Missouri  basin  is  comprised  of  relatively  flat  plains  while  most  of  the  Ohio  basin  is  more 
mountainous.  The  vegetation  types  would  also  be  different  between  the  two  basins  The 
model  runoff  coefficients  suggest  that  the  evaporation  of  water  over  the  Missouri  basin  is 
too  low  and  the  evaporation  over  the  Ohio  River  is  too  high.  The  difference  in  the 
observed  and  model  runoff  coefficients  for  the  Missouri  and  Ohio  rivers  balance  each 
other  out  closer  to  the  mouth  of  the  Mississippi  basin  where  the  observed  and  model 
runoff  coefficients  for  the  Mississippi  River  at  Vicksburg  are  .22  and  .23,  respectively 

Table  4.2  shows  the  model-generated  snow  mass  for  the  various  basins  averaged 
for  the  fall,  winter,  and  spring  seasons.  Most  snow  falls  in  the  Missouri  basin,  which 
produces  a  river  runoff  maximum  in  April  and  May.  However,  the  observed  river  runoff  is 
much  lower  and  shows  little  monthly  variation  suggesting  that  the  model  generates  too 
much  snow  and  snow  melt  runoff  for  the  Missouri  River.  The  model  generates  very  low 
snow  amounts  for  the  Ohio  River  suggesting  model-generated  river  runoff  there  is 
dependent  on  precipitation  only.  Observed  river  runoff  in  the  Ohio  River  shows  a 
maximum  in  March  despite  uniform  precipitation  throughout  the  year  This  suggests  that 
snow  melt  is  contributing  to  the  observed  runoff  maximum.  The  model  generates  too  little 
snow  in  the  Ohio  basin  causing  less  runoff  than  observed  and  generates  too  much  snow  in 
the  Missouri  basin  causing  more  runoff  than  observed. 
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Table  4.2.  Seasonal  model-generated  snow  mass  (Km-*) 
averaged  over  fall,  winter,  and  spring  for  the  Mississippi 
River  and  its  tributaries  for  the  present  climate _ 


Average  Seasonal  Snow  Mass 
Present  Climate 

River  Basin 

Fall 

Winter 

Spring 

Illinois 

0.0 

0.4 

0.0 

Missouri 

1.4 

39.0 

86 

Ohio 

0.0 

0.7 

0.0 

Arkansas 

0.2 

1.7 

0.1 

Upper  Mississippi 

0.4 

6.6 

0.6 

Central  Mississippi 

1.8 

46.2 

9.3 

Lower  Mississippi 

2.1 

50.2 

9.7 

One  problem  can  be  the  interpolation  of  precipitation  and  snow  mass  from  a  4°  X 
5°  grid  box  to  four  2°  X  2.5°  grid  boxes.  It  is  assumed  that  precipitation  occurs  uniformly 
over  the  entire  grid  box.  The  problem  with  the  Missouri  River  basin  is  with  the  northwest 
grid  boxes  bordering  the  mountains.  The  most  western  2°  X  2.5°  grid  boxes  were 
interpolated  from  a  4°  X  5°  grid  box  that  covered  part  of  a  mountainous  region  Only  half 
of  the  4°  X  5°  grid  box  was  within  the  Missouri  River  basin  Because  precipitation  and 
snow  fall  is  assumed  uniform  over  the  entire  4°  X  5°  grid  box  it  was  the  same  for  each  2° 
X  2.5°  grid  box 

Another  related  problem  is  that  grid  boxes  may  contain  several  rivers  which  belong 
to  different  basins.  Water  from  the  entire  grid  box  is  assumed  to  flow  in  one  direction, 
thus,  river  runoff  may  enter  the  wrong  river  basin.  This  is  also  happening  in  the  Missouri 
River  basin  in  which  the  western  most  grid  box  contains  precipitation  and  snow 
interpolated  from  a  4  X  5  grid  box  encompassing  both  the  Missouri  and  Snake  rivers. 
Since  grid  box  runoff  must  be  assigned  to  a  particular  basin,  some  runoff  intended  to  flow 
into  the  Snake  River  is  flowing  into  the  Missouri  River. 

The  goal  of  this  chapter  was  to  further  analyze  a  large  river  basin  by  dividing  it 
into  its  tributaries.  This  could  be  done  for  any  of  the  rivers  discussed  in  chapter  3.  The 
Mississippi  River  was  chosen  because  annual  precipitation  was  simulated  well  while  the 
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model-generated  monthly  precipitation  differed  in  the  spring  and  fall  and  good  observed 
data  was  available.  The  model  generates  too  much  precipitation  in  the  western  Mississippi 
and  too  little  in  the  eastern  Mississippi  causing  differences  in  the  simulations  of  monthly 
and  annual  river  runoff.  Care  must  be  given  to  areas  of  the  grid  which  border 
mountainous  regions.  Because  precipitation  and  snow  amounts  are  assumed  uniform 
throughout  a  grid  box,  interpolation  into  smaller  grid  boxes  can  cause  errors  in 
precipitation  amounts  and  eventually  in  the  routing  of  the  runoff.  For  the  Missouri  and 
Ohio  rivers,  model-generated  snow  seems  to  be  reversed,  with  the  model  generating  too 
much  snow  for  the  Missouri  and  too  little  snow  for  the  Ohio  River  However,  this  study 
did  not  utilize  observed  snow  data  for  the  Missouri  and  Ohio  basins  causing  this 
conclusion  Finer  spatial  and  temporal  resolution  should  improve  the  model's  ability  to 
simulate  river  runoff  and  precipitation  for  the  world's  river  basins 
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Chapter  5.  Changes  in  monthly  runoff  in  a  doubled  CO2  climate 

Future  changes  in  river  runoff  will  impact  many  areas,  including  agriculture,  water 
resources,  and  land  use.  In  this  chapter  potential  changes  in  river  runoff  are  examined  for 
a  doubled  CO2  climate.  Miller  and  Russell  [1992]  examined  mean  annual  river  runoff  for 
the  present  climate  and  for  a  doubled  CO2  climate  and  found  increased  river  runoff  for  25 
of  the  world's  33  major  rivers.  The  largest  increases  were  found  in  the  higher  latitudes, 
where  substantial  shifts  in  runoff  patterns  may  occur  because  temperature  changes  can 
effect  the  ratio  of  rain  to  snow  Less  snow  during  the  winter  could  have  a  major  impact 
on  water  resources  during  spring  and  summer.  This  chapter  extends  the  work  of  Miller 
and  Russell  [1992]  to  examine  seasonal  changes  in  river  runoff  with  emphasis  on  high 
latitudes  where  snow  melt  is  an  important  component  of  river  runoff  It  is  important  to 
examine  monthly  runoff  to  determine  whether  changes  are  uniform  throughout  the  year  or 
are  caused  by  seasonal  variations  in  precipitation  and  river  runoff 

The  off-line  method  discussed  in  chapter  1  can  be  used  with  temperatures  and 
precipitation  from  GCMs  or  other  sources.  Gleick  [1987]  and  Leftenmaier  and  Gan 
[1990]  have  used  precipitation  and  temperature  from  GCMs.  Gleick  [1987]  examined  the 
Sacramento  River  basin  in  California  for  18  widely  varying  climate  changes,  ten  from 
hypothetical  temperature  and  precipitation  changes  and  eight  from  precipitation  and 
temperature  changes  generated  by  GCMs.  He  found  that  climate  change  caused  increased 
runoff  in  the  winter  and  decreased  runoff  in  the  summer.  The  principal  physical 
mechanism  concluded  by  Gleick  [1987]  was  a  decrease  snow.  Leftenmaier  and  Gan 
[1990]  found  similar  results  for  several  northern  California  river  basins. 

Flaschko  el  al.  [1987]  used  the  off-line  method  and  examined  changes  of  river 
runoff  in  the  Great  Basin  Region  of  the  western  United  States  by  applying  water  balance 
models  to  four  watersheds  in  Nevada  and  Utah  They  modeled  the  effects  of  four  climatic 
change  scenarios  suggested  by  Stockton  and  Boggess  [1979],  1)  a  2°C  increase  in 
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temperature  and  a  10%  decrease  in  precipitation  (warm/dry),  2)  a  2°C  increase  in 
temperature  and  a  25%  decrease  in  precipitation  (warm/very  dry),  3)  a  2°C  decrease  in 
temperature  and  a  10%  increase  in  precipitation  (cool/wet),  and  4)  a  2°C  decrease  in 
temperature  and  a  25%  increase  in  precipitation  (cool/very  wet),  on  the  Carson,  Martin, 
Bear,  and  Sevier  basins,  with  emphasis  on  the  warmer/dry  case,  the  more  likely  scenario. 
They  found  that  runoff  decreased  by  1 7%  to  28%  from  the  present  mean  for  the  four  river 
basins.  Roos  [1989]  and  Ayers  et  al.  [1990]  used  the  same  type  of  water  balance  model  to 
examine  changes  in  runoff  for  the  Sacramento  River,  San  Joaquin  River,  and  rivers  in  the 
Tulare  Lake  system  of  California  and  the  Delaware  River,  respectively,  with  emphasis  on 
spring  snow  melt.  In  both  instances,  precipitation  patterns  were  assumed  to  remain  the 
same  in  the  doubled  CO2  climate,  however,  because  of  the  wanner  temperatures,  the 
majority  of  the  winter  precipitation  was  in  the  form  of  rain  with  much  less  snow  in  the 
mountain  regions.  Thus,  spring  runoff  was  reduced  impacting  both  agriculture  and  water 
resources. 

This  chapter  will  focus  on  river  runoff  using  grid  box  runoff  produced  directly 
from  a  GCM.  Three-year  simulations  for  the  present  climate  and  the  doubled  CO2  climate 
were  performed  at  4°  X  5°  resolution  using  the  GISS  model.  The  same  simulations  were 
used  in  Miller  and  Russell  [1992]  for  examining  changes  in  mean  annual  runoff.  As 
discussed  in  Miller  and  Russell  [1992],  this  relatively  short  simulation  introduces  the 
potential  for  temporal  sampling  errors.  Such  errors  will  be  greater  in  monthly  computed 
runoff.  The  annual  cycle  of  sea  surface  temperatures  (SST)  and  ice  distribution  were 
specified  by  interpolating  from  equilibrium  simulations  at  8°  X  10°  resolution,  in  which 
SST  and  ice  distribution  were  predicted  [Hansen  et  al. ,  1984].  The  model  simulations  are 
A5 1  and  947  as  shown  in  Table  3.1. 

The  model-generated  annual  river  runoff  and  precipitation  for  the  present  climate 
and  doubled  CO2  climate  are  shown  in  Table  5.1  and  Figures  5. 1-5. 2  for  30  river  basins 
These  results  are  not  identical  to  those  of  Miller  and  Russell  [  1 992]  because  the  grid 
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boxes  assigned  to  specific  rivers  have  been  modified  somewhat  and  observed  runoff  for 
this  study  is  from  UNESCO  [1969,  1974,  1985]  and  not  from  Millman  and  Meade  [1983] 
Drainage  areas  for  the  model  are  compared  with  the  observed  areas  in  Table  5  2 
Observed  river  runoff  and  areas  are  from  UNESCO  [1969,  1974,  1985]  while  observed 
precipitation  is  from  Legates  and  Willmott  [1990],  Observed  areas  for  runoff  may  differ 
from  the  model-generated  areas  because  the  observed  stations  are  not  always  at  the  river 
mouth.  Hence,  in  most  cases  the  model-generated  runoff  is  based  on  a  larger  drainage 
basin  area  than  the  observed 

The  world's  river's  were  divided  into  three  categories,  those  with  less  than  60  cm 
of  annual  precipitation,  those  with  between  60  and  120  cm,  and  those  with  more  than  120 
cm  yr1  which  are  referred  to  as  dry,  moderately  wet,  and  wet  [Kuhl  and  Miller  1992] 
Dry  basins  are  separated  into  dry  and  high  latitude  dry  (north  of  45°N)  basins  When 
averaged  over  all  30  rivers,  the  model-generated  annual  precipitation  and  runoff  are  25% 
and  13%  too  high,  respectively,  for  the  present  climate.  Model-generated  annual 
precipitation  is  high  for  all  four  categories  with  the  greatest  discrepancy  in  the  dry  river 
basins  where  the  model-generated  precipitation  is  more  than  50%  too  high  This  is  also 
true  for  the  annual  model -generated  river  runoff,  where  the  model  does  poorest  in  the  dry 
river  basins.  Rind  et  al.  [1990]  suggests  that  (3  in  Eq.  2.7  is  not  calculated  accurately  in 
the  GISS  model.  This  could  lead  to  incorrect  evaporation  and  incorrect  runoff.  The 
model  also  does  not  allow  for  runoff  to  evaporate  as  it  moves  downstream  This  effect 
would  be  most  pronounced  in  the  dry  basins. 

Figures  5.3a-5.3d  show  the  seasonal  changes  (Km3)  in  the  model-generated 
precipitation  and  runoff  for  the  doubled  CO2  climate  The  greatest  changes  in 
precipitation  for  the  wet  climates  (Fig  5-3a)  occurs  between  September  and  February 
The  Congo  and  to  a  lesser  extent  the  Amazon  are  the  only  rivers  in  which  the  precipitation 
increases  are  fairly  uniform  for  all  seasons.  The  smallest  changes  occur  between  March 
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and  May.  The  Magdalena  is  the  only  river  to  show  decreases  in  both  nrnoff  and 
precipitation. 

The  Niger  River  (Fig  5-3  c)  is  particularly  interesting  because  an  increase  in 
precipitation  from  December  through  May  is  accompanied  by  a  decrease  in  runoff  for  the 
same  period.  The  decreased  runoff  is  further  enhanced  in  the  latter  part  of  the  year  when  a 
decrease  in  precipitation  also  occurs.  The  decrease  in  runoff  is  due  to  increased 
evapotranspiration. 

High  latitude  rivers  (Fig  5-3  d)  show  the  greatest  change  in  runoff  to  occur  in  the 
spring  months  (MAM),  due  tc  a  combination  of  snow  melt  and  increase  precipitation  The 
confidence  in  these  results  is  reduced  somewhat  because  the  model-generated  river  runoff 
and  precipitation  are  generally  too  high  in  the  present  climate  and  also  because  the 
simulations  are  for  only  three  years.  The  short  simulations  are  likely  to  be  more  of  a 
factor  for  the  Magdalena  and  Sao  Francisco  rivers  because  the  drainage  areas  are  small.  A 
more  detailed  examination  of  severat  river  basins  follows. 

Figures  5.4-5.16  show  monthly  precipitation  and  river  runoff  for  the  present  and 
doubled  COj  climates.  The  world's  two  largest  rivers,  the  Amazon  and  Congo,  (Fig  5  4) 
are  classified  as  wet  river  basins.  Results  for  these  two  rivers  are  not  too  reliable  because 
river  runoff  is  poorly  simulated  for  the  present  climate  The  Amazon  has  increased 
precipitation  between  November  and  February,  the  southern  hemisphere  summer,  and 
increased  runoff  from  November  through  April  As  discussed  in  chapter  3,  model¬ 
generated  river  runoff  using  the  new  routing  scheme,  will  peak  after  the  precipitation 
maximum  by  approximately  1  to  2  months.  The  routing  scheme  allows  for  movement  of 
runoff  from  one  grid  box  to  another  until  it  reaches  the  mouth  of  the  river  basin  The 
Congo  river  has  increased  precipitation  and  runoff  throughout  the  entire  year,  with  a  slight 
maximum  in  the  spring  and  summer  months. 

The  Sao  Francisco  river  (Fig.  5  5)  is  one  of  the  rivers  in  which  annual  precipitation 
increases  slightly  (6%)  while  annual  runoff  decreases  slightly  (3%)  Monthly  precipitation 
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and  runoff  show  little  change  throughout  the  year.  The  Mekong  River  is  shown  in  Fig  5  6 
Besides  precipitation  and  runoff,  average  monthly  snow  mass  is  shown  for  the  snow 
season.  The  snow  decrease  is  primarily  due  to  the  warmer  temperatures  associated  with 
the  doubled  CO2  climate.  Since  we  have  not  compared  model-generated  snow  fall  with 
observed  snow  fall,  this  result  is  tentative.  Because  snow  decreases,  one  would  expect  a 
decrease  in  snow  melt  runoff  during  the  late  spring  and  summer  months.  However,  the 
Mekong  river  runoff  actually  increases  during  these  months  because  of  the  increased 
spring  and  summer  rainfall. 

Figures  5. 7-5.9  show  six  of  the  basins  classified  as  moderately  wet.  Figure  5.5 
shows  the  Mississippi  and  Yangtze  rivers.  Changes  in  snow  mass  can  significantly  affect 
river  runoff  Although  snow  mass  decreases  for  both  rivers  during  the  winter  months, 
spring  runoff  shows  little  change  because  of  the  increased  rainfall  during  the  spring  and 
summer 

The  Niger  (Fig  5.8)  and  the  Nile  (Fig.  5.9)  rivers  are  two  moderately  wet  rivers  in 
which  there  is  little  change  in  annual  precipitation.  However,  a  decrease  in  precipitation 
occurs  from  June  to  October  and  an  increase  in  precipitation  occurs  from  November  to 
February  for  the  Nile  River  and  November  through  April  for  the  Niger  River  The 
precipitation  decrease  between  June  and  October  is  approximately  equal  for  both  the  Nile 
and  Niger  rivers,  however,  the  decrease  in  river  runoff  is  much  greater  for  the  Niger  River. 
The  decreased  river  runoff  must  be  due  to  increased  evapotranspi ration.  Observed 
temperatures  average  26°  -  27°C  while  precipitation  varies  significantly  over  the  Niger 
basin  Temperature  increases  in  response  to  global  warming  may  only  enhance 
evaporation  in  the  drier  regions  of  the  Niger  River.  Unfortunately  the  observed  station  for 
the  Niger  River  is  far  away  from  the  mouth  of  the  river.  Hence,  one  cannot  be  certain 
about  the  model's  ability  to  simulate  the  runoff  for  the  present  climate. 

The  model -generated  river  runoff  for  the  low  and  mid-latitude  dry  river  basins  is 
poorly  simulated.  Hence,  the  confidence  in  the  runoff  results  for  these  basins  is  much 
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lower  than  for  the  other  cases.  The  model-generated  runoff  for  the  Colorado  River  (Fig 
5.10b)  is  too  high  when  compared  to  observed  runoff  and  precipitation  is  poorly 
simulated  The  major  change  in  the  precipitation  for  the  Yellow  River  (Fig  5  10a)  is  an 
enhancement  of  the  summer  monsoon  in  the  doubled  CO2  climate  River  runoff  is 
maximum  in  October,  approximately  one  month  after  the  monsoon  season. 

River  runoff  at  the  higher  latitudes  is  more  interesting  to  analyze  because  of  the 
contribution  of  snow  melt.  Because  observed  snow  data  have  not  been  obtained  to 
compare  with  the  model-generated  snow,  the  ability  of  the  model  to  simulate  the  actual 
snow  mass  is  not  known.  However,  changes  in  the  total  snow  mass  in  a  doubled  C02 
climate  can  be  used  to  further  analyze  the  effects  of  snow  melt  on  river  runoff  in  the  higher 
latitude  rivers.  Table  5.3  shows  the  seasonal  average  snow  mass  for  both  the  present 
climate  and  the  doubled  C02  climate.  For  all  rivers,  except  the  Amur,  which  shows  little 
change,  fall  snow  mass,  in  a  doubled  C02  climate,  is  reduced  by  an  average  of  30%. 

There  is  an  interesting  difference  between  northeast  Asia  and  other  high  latitude 
rivers.  In  northeast  Asia  the  net  snow  mass  actually  increases  for  the  Yenesei,  Lena, 
Kolyma,  Indigirka,  and  Amur  rivers.  For  all  other  high  latitude  rivers,  the  average  winter 
and  spring  snow  mass  decrease.  The  snow  decrease  is  caused  primarily  by  an  increase  in 
temperatures  for  a  doubled  C02  climate.  However,  because  the  average  temperatures  in 
the  northeastern  region  of  Asia  are  well  below  freezing,  temperature  changes  of  2-7 
degrees  Celsius  will  not  cause  changes  in  the  rain  to  snow  ratio,  thus,  an  increase  in 
precipitation  in  a  doubled  C02  climate  will  generate  more  snow  mass  for  the  river  basins 
in  that  region.  Despite  differences  in  amounts  of  snow  mass,  high  latitude  rivers  all  show 
one  thing  in  common,  a  faster  snow  melt  during  spring  in  a  doubled  C02  climate.  This 
faster  snow  melt  causes  the  river  runoff  maximum  to  occur  earlier  for  ail  high  latitude  river 
basins 

The  two  high  latitude  rivers  in  North  America  are  the  Yukon  and  the  Mackenzie 
(Fig.  5.12).  For  both  rivers  there  is  a  decrease  in  the  monthly  snow  mass  throughout  the 
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snow  season,  and  the  melting  season  begins  approximately  one  month  earlier  in  the 
doubled  CO2  climate  This  causes  a  shift  in  the  timing  of  the  spring  river  runoff  peak 
However,  the  magnitude  of  the  maximum  river  runoff  is  the  same  for  the  Yukon  and 
increases  slightly  for  the  Mackenzie  due  to  increased  spring  precipitation  in  the  doubled 
CO2  climate.  The  shift  in  the  runoff  peak  also  occurs  in  the  Sevemay  Dvina  and  Ob  rivers 
located  in  north  central  Asia  (Fig.  5.13).  The  changes  in  timing  of  the  maximum  river 
runoff  could  be  studied  better  if  daily  data  were  available. 

Figures  5.14-5.16  show  the  Yenesei,  Lena,  Kolyma,  Indigirka,  and  Amur  rivers,  all 
of  which  have  an  increase  in  average  monthly  snow  mass  during  the  winter  and  spring 
months.  However,  despite  the  increase  in  snow  mass,  the  snow  melt  ends  at  about  the 
same  time  for  both  the  present  and  doubled  CO2  climates.  The  increase  in  snow  mass  is 
associated  with  increased  winter  precipitation,  however,  unlike  the  other  high  latitude 
rivers,  temperature  increases  are  not  large  enough  to  push  the  temperature  over  the 
freezing  point.  Climatology  for  the  present  climate  from  Critchfield  [1974]  shows  high 
pressure  is  centered  over  the  Tibetan  Plateau,  with  the  average  surface  flow  to  be 
southerly  to  southwesterly  for  all  river  basins  west  of  the  Yenesei  for  January  The 
Yenesei  river  basin  and  river  basins  to  the  east  are  under  an  average  westerly  to 
northwesterly  surface  flow  for  January.  Observed  January  temperatures  average  between 
-50°C  and  -30°C  for  the  regions  east  of  the  Yenesei  river  basin  and  between  -25°C  and  - 
10°C  for  the  region  west  of  the  Yenesei  river  basin.  In  the  doubled  C02  climate,  the 
temperature  is  projected  to  rise  between  2°C  and  7°C  for  the  high  latitude  rivers.  The 
temperature  increase  is  not  affecting  the  snow  to  rain  ratio  in  river  basins  in  which  the 
average  temperatures  are  well  below  freezing  This  can  explain  the  differences  in  snowfall 
between  the  high  latitude  river  basins. 

The  major  change  in  the  high  latitude  rivers  is  that  the  spring  runoff  peak  occurs 
one  month  earlier  in  a  doubled  CO2  climate  The  shift  is  primarily  due  to  the  earlier 
melting  of  snow  and  ice  Despite  an  increase  in  precipitation  for  all  high  latitude  rivers. 
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the  magnitude  of  the  runoff  peak  is  greater  in  the  doubled  CO2  climate  for  only  the 
Yenesei,  Lena,  Kolyma,  Indigirka,  and  Amur  rivers.  This  is  caused  by  the  combination  of 
increased  precipitation  and  snow. 

One  advantage  of  the  new  routing  scheme  of  Miller  et  al.  [1992]  is  that  larger 
river  basins  can  be  divided  into  smaller  basins.  The  Mississippi  River  and  its  tributaries 
were  examined  for  a  doubled  CO2  climate.  Figures  5.17-  5-20  show  present  and  doubled 
CO2  climates  for  precipitation,  runoff,  and  snow  mass.  Annual  runoff  and  precipitation 
increase  16%  and  21%,  respectively,  and  snow  mass  decreased  for  the  Mississippi  River 
basin.  The  most  significant  changes  occurred  in  the  Missouri  basin  in  which  runoff 
maxima  occur  in  March  and  May.  The  combination  of  decreases  in  snow  mass  and  April 
precipitation  cause  this  double  maxima  to  occur.  The  confidence  in  these  results  is  not  too 
high  because  the  major  tributaries,  the  Missouri  and  Ohio,  are  poorly  simulated  by  the 
model.  For  the  present  climate  the  model  generates  too  much  precipitation  and  runoff  for 
the  Missouri  and  too  little  for  the  Ohio 

Annual  and  monthly  changes  in  precipitation  and  river  runoff  in  a  doubled  C02 
climate  were  examined  for  several  river  basins  worldwide.  It  is  important  to  examine 
monthly  changes  to  determine  whether  any  of  the  changes  are  uniform  throughout  the 
year  Most  river  basins  have  increased  precipitation  and  river  runoff  in  a  doubled  C02 
climate.  The  one  exception  is  the  Niger  river  in  which  their  is  little  change  in  precipitation 
and  a  reduction  in  river  runoff  of  36%,  caused  by  increased  evapotranspiration  The  most 
interesting  changes  in  runoff  occur  in  the  higher  latitudes  where  snow  melt  is  part  of  the 
total  river  runoff.  The  spring  peak  is  earlier  and  the  magnitude  is  higher  for  the 
northeastern  rivers  in  Asia  where  snow  mass  increases  in  a  doubled  C02  climate.  The 
increase  in  snow  mass  is  unique  to  only  the  rivers  in  northeastern  Asia 
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Table  5.1.  Annual  observed  and  model-generated  runoff  (Km-Vyr)  and  precipitation  (Km^/yr)  for  the 
world’s  major  rivers  for  the  present  climate  (1XC02)  and  the  doubled  CO2  climate  (2XC02).  Change 
from  IXCO2  to  2XCO2  is  given  by  the  %  change  column. 


Climate  Type 
River 

Precipitation 

OBS  1XCQ2  2XC02 

OBS 

Runoff 

1XC02 

2XC02 

%  Change 
Precip  Runoff 

Wet 

Amazon 

13942 

14224 

15538 

4886 

3079 

3341 

9% 

9% 

Congo 

5516 

8470 

9766 

1398 

2095 

2632 

15% 

26% 

Orinoco 

2193 

2005 

2422 

794 

466 

682 

21% 

46% 

Mekong 

1477 

2138 

2317 

449 

711 

815 

8% 

15% 

Magdalena 

758 

1056 

924 

213 

303 

217 

-13% 

-28% 

Sao  Francisco 

799 

1273 

1344 

83 

274 

265 

6% 

-3% 

Average 

4114 

4861 

5385 

1304 

1155 

1325 

11% 

15% 

Moderately  Wet 

Yangtze 

1985 

3465 

4126 

792 

1499 

1947 

19% 

30% 

Mississippi 

3039 

2981 

3559 

498 

660 

766 

19% 

16% 

LaPlata  (Parana) 

3662 

3164 

4064 

— 

(475) 

(723) 

28% 

52% 

St  Lawrence 

1070 

1102 

1256 

214 

321 

350 

14% 

9% 

Danube 

1192 

1683 

1984 

— 

(516) 

(585) 

18% 

13% 

Columbia 

507 

774 

827 

172 

249 

265 

7% 

6% 

Zambesi 

1218 

1356 

1475 

105 

197 

213 

9% 

8% 

Fraser 

179 

307 

329 

87 

134 

146 

7% 

9% 

Nile 

1960 

3659 

3690 

— 

(657) 

(633) 

1% 

-4% 

Niger 

1618 

2116 

2141 

— 

(332) 

(211) 

1% 

-36% 

Average 

1643 

2061 

2345 

311 

510 

615 

14% 

16% 

Dry 

Indus 

392 

592 

489 

76 

204 

124 

-17% 

-39% 

Tigris-Euphrates 

497 

497 

481 

— 

(76) 

(77) 

-3% 

1% 

Yellow 

541 

1454 

1642 

— 

(563) 

(710) 

13% 

26% 

Colorado 

205 

339 

385 

12 

69 

91 

14% 

32% 

Murray 

528 

585 

702 

8 

97 

122 

20% 

26% 

Average 

433 

693 

740 

32 

202 

225 

7% 

11% 

Hieh  Latitude  Dry 

Yenesei 

1126 

1371 

1770 

558 

462 

651 

29% 

41% 

Lena 

926 

1486 

1712 

516 

535 

681 

15% 

27% 

Ob 

1354 

1066 

1564 

388 

464 

632 

47% 

36% 

Amur 

1120 

1394 

1438 

309 

273 

298 

3% 

9% 

Mackenzie 

607 

1242 

1504 

264 

644 

756 

21% 

17% 

Yukon 

333 

715 

824 

197 

511 

630 

15% 

23% 

Sevemay  Dvina 

206 

241 

322 

107 

132 

165 

34% 

25% 

Kolyma 

258 

506 

692 

71 

367 

526 

37% 

43% 

Indigirka 

100 

228 

275 

49 

134 

169 

21% 

26% 

Average 

670 

917 

1122 

273 

391 

501 

22% 

28% 

Total  Average 

1644 

2050 

2319 

510 

578 

687 

13% 

18% 

Observed  precipitation  is  from  Legates  and  WiUmott  [1990]  Observed  runoff  ts  from  UNESCO  [1969,  1974, 
1985].  ( — )  Indicates  river  basins  in  which  observed  runoff  station  is  not  near  the  mouth  of  the  nver.  Values  in 
(  )  are  not  included  in  climate  type  and  total  runoff  averages. 
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Table  5.2.  Model  and  observed  drainage  basin  areas  (Km2)  for  the  world’s  major  rivers  The  ratio  is 
defined  as  the  model  to  observed  ratio. 


Climate  Type 
River 

Model 

Area 

Observed 

Ratio 

Wet 

Amazon 

6500000 

4640285 

1.40 

Congo 

3510000 

3475000 

1.01 

Orinoco 

950000 

850000 

1.12 

Mekong 

860000 

646000 

1.33 

Magdalena 

360000 

257438 

1.40 

Sao  Francisco 

650000 

622600 

1.04 

Average 

2138333 

1748554 

1.22 

Moderatelv  Wet 
Yangtze 

1810000 

1705383 

1  06 

Mississippi 

3510000 

2964300 

1.18 

LaPlata  (Parana) 

2910000 

975375 

2.98 

St  Lawrence 

1140000 

764600 

1.49 

Danube 

1530000 

807000 

1.90 

Columbia 

710000 

614000 

1.16 

Zambesi 

1190000 

940000 

1.27 

Fraser 

230000 

217000 

1.06 

Nile 

2760000 

— 

— 

Niger 

1490000 

— 

— 

Average 

1728000 

1123457 

1.54 

Dry 

Indus 

830000 

832418 

1  00 

Tigris-Euphrates 

1150000 

408100 

2.82 

Yellow 

1060000 

688421 

1.54 

Colorado 

640000 

629100 

1.02 

Murray 

1040000 

991000 

1.05 

Average 

944000 

709808 

1.33 

High  Latitude  Drv 

Yenesei 

2700000 

2440000 

111 

Lena 

2400000 

2430000 

0.99 

Ob 

2630000 

2430000 

1.08 

Amur 

1870000 

1730000 

1.08 

Mackenzie 

1570000 

1570000 

1.00 

Yukon 

770000 

767000 

1.00 

Sevemay  Dvina 

340000 

348000 

0.98 

Kolyma 

710000 

361000 

1.97 

Indigirka 

350000 

305000 

1.15 

Average 

1482222 

1375667 

1.08 

Total  Average _ 

1605667 

1264608 

1.27 

Observed  area  from  ItNESCO  [1969,  1974,  1985] 
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Table  5.3.  Seasonal  model-generated  snow  mass  (Km-*)  averaged  over  fall,  winter,  and  spring  for  the 
world's  major  rivers  for  the  present  climate  (1XC02)  and  the  doubled  CO2  climate  (2XC02).  Change 
from  IXCO2  to  2XCO2  is  given  by  the  %  change  column 


Climate  Type 
River 

Fall 

1XC02 

Fail 

2XC02 

Average  Seasonal  Snow  Mass 

Fail  Winter  Winter  Winter 

•4  Change  1XC02  2XC02  %  Change 

Spring 

1XC02 

Spring 

2XC02 

Spring 

%  Change 

Wet 

Mekong 

6.7 

1.7 

-75*/. 

21.3 

18.3 

-14% 

17.3 

11.0 

-36% 

Average 

6.7 

1.7 

-75% 

21.3 

18.3 

-14% 

17.3 

11.0 

-36% 

Moderately  Wet 
Yangtze 

12.0 

3.3 

-73% 

54.3 

28.3 

-48% 

31.7 

16.7 

-47% 

Mississippi 

2.0 

0.7 

-65% 

49.7 

19.7 

-60% 

97 

1.3 

-87% 

St  Lawrence 

3.3 

1.0 

-70% 

58.7 

24  3 

-59% 

36.3 

10.3 

-72% 

Danube 

0.3 

0.0 

--- 

17.0 

1.7 

-90% 

2.3 

0.0 

— 

Columbia 

0.3 

0.0 

... 

8.3 

4.3 

-48% 

1.0 

0.0 

— 

Fraser 

0.3 

0.0 

... 

8.3 

2.0 

-76% 

4.7 

0.0 

... 

Average 

3.0 

0.8 

-73% 

32.7 

13.4 

-59% 

14  3 

4.7 

-67% 

Dry 

Indus 

00 

0.0 

... 

14.7 

00 

--- 

5.3 

0.0 

Yellow 

1.7 

0.0 

... 

7.0 

3.7 

-47% 

2.3 

0.3 

-87% 

Colorado 

0.0 

0.0 

... 

4.0 

0.7 

-83% 

0.7 

0.0 

... 

Average 

0.6 

0.0 

... 

8.6 

1.5 

-83% 

2.8 

0.1 

-96% 

High  Latitude  Dry 

Yenesei 

34.0 

25.3 

-26% 

221  0 

247.0 

7% 

160.7 

If  7.3 

17% 

Lena 

44.3 

35.3 

-20% 

203.0 

225.0 

11% 

187  7 

205,3 

9% 

Ob 

22.0 

14.0 

-36% 

205.0 

171.3 

-16% 

118  7 

88.0 

-26% 

Amur 

6.0 

6.3 

5% 

62.3 

68.7 

10% 

35.7 

39.0 

9% 

Mackenzie 

31.3 

21.0 

-33% 

216.0 

179.3 

-17% 

212.0 

134.0 

-37% 

Yukon 

22.7 

16.3 

-28% 

129.3 

103.0 

20% 

166.7 

112.3 

-33% 

Sevemay  Dvina 

5.3 

1.7 

68% 

45.3 

29.3 

-35% 

300 

12  0 

-60% 

Kolyma 

26.3 

20.3 

-23% 

108.0 

125  3 

16% 

148,7 

170.C 

14% 

Indigirka 

140 

10.7 

-24% 

43  7 

46  7 

7% 

60.0 

62.3 

4% 

Average 

22.9 

16.8 

-27% 

137.1 

131.7 

-4% 

124.5 

112.2 

-10% 

Total  Average 

12.2 

8.3 

-32% 

77.7 

67.8 

-13% 

64.8 

55.3 

-15% 

Seasons  are  defined  as  the  following:  Fall  (September,  October,  and  November),  Winter  (December,  January  , 
and  February),  and  Spring  (March,  April,  and  May).  The  amounts  represent  the  average  snow  mass  over  the 
three  month  period. 


runoff  seasonal  changes  for  the  wet  river  basins  in  the  doubled 
Aug  (JJA),  and  Sep,  Oct,  Nov  (SON). 


Precipitation  CJ Runoff  I  j  L_j  Precipitation 


Figure  5.3c.  Model-generated  precipitation  and  runoff  seasonal  changes  for  the  dry  river  basins  in  the  doubled  CO2  climate  for  Dec, 
Jan,  Feb  (DJF),  Mar,  Apr,  May  (MAM),  Jun,  Jul,  Aug  (JJA),  ar  ’  jep,  Oct,  Nov  (SON). 


Precipitation 


Figure  5  4  Model-generated  monthly  precipitation  and  runoff  for  the  present  climate  [M(])]  and  the  doubled  CO2  climate  [M(2)]  and 
observed  [O]  precipitation  from  Legates  and  Willmott  [1990]  and  runoff  from  UNESCO  [1969,  1974,  1985]  for  the  (a)  Amazon  and 
(b)  Congo  rivers 


MEKONG  MEKONG  MEKONG 
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Figure  5.6  Model-generated  monthly  precipitation,  runoff,  and  averaged  snow  mass  during  the  snow  season  for  the  present  climate 
[M(l)]  and  the  doubled  CO2  climate  [M(2)]  and  observed  [O]  precipitation  from  Legates  and  Willmott  [1990]  and  runoff  from 
UNESCO  [1969,  1974,  1985]  for  the  Mekong  River. 


NIGER  NIGER 
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Figure  5.8  Model-generated  monthly  precipitation  and  runoff  for  the  present  climate  [M(l)]  and  the  doubled  CO2  climate  [M(2)]  and 
observed  [O]  precipitation  from  Legates  and  Willmott  [1990]  and  runoff  from  UNESCO  [1969,  1974,  1985]  for  the  (a)  Niger  and  (b) 
LaPlata  (Parana)  rivers.  Observed  runoff  stations  for  Niger  and  LaPlata  are  not  near  the  mouth  of  the  river 


ZAMBESI  ZAMBESI 
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Figure  5  9  Model-generated  monthly  precipitation  and  runoff  for  the  present  climate  [M(l)]  and  the  doubled  C02  climate  [M(2)]  and 
observed  [O]  precipitation  from  Legates  and  Willmott  [1990]  and  runoff  from  UNESCO  [1969,  1974,  1985]  for  the  (a)  Zambesi  and 
(b)  Nile  rivers.  Observed  runoff  station  for  the  Nile  is  not  near  the  mouth  of  the  river. 


TIGRIS  EUPHRATES  TIGRIS  EUPHRATES 
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Figure  5.11.  Model-generated  monthly  precipitation  and  runoff  for  the  present  climate  [M(l)j  and  the  doubled  COj  climate  [M(2)]  and 
observed  [O]  precipitation  from  Legates  and  Willmott  [1990]  and  runoff  from  UNESCO  [1969,  1974,  1985]  for  the  (a)  Tigris- 
Euphrates  and  (b)  Murray  rivers  Observed  runoff  station  for  the  Tigris-Euphrates  is  not  near  the  mouth  of  the  river. 


Precipitatin' 
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Figure  5. 12.  Model-generated  monthly  precipitation,  runoff,  and  averaged  snow  mass  during  the  snow  season  for  the  present  climate 
[M(l)]  and  the  doubled  C02  climate  [M(2)]  and  observed  [O]  precipitation  from  Legates  and  WUlmott  [1990]  and  runoff  from 
UNESCO  [1969.  1974.  19851  for  the  (a)  Yukon  and  (b)  Mackenzie  rivers 
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Precipitation 


Figure  5  15  Model-generated  monthly  precipitation,  runoff,  and  averaged  snow  mass  during  the  snow  season  for  the  present  climate 
[M(l)]  and  the  doubled  CO2  climate  [M(2)]  and  observed  [O]  precipitation  from  Legates  and  Willmott  [1990]  and  runoff  from 
UNESCO  [1969,  1974,  1985]  for  the  (a)  Kolyma  and  (b)  Indigirka  rivers. 


UPPER  MISSISSIPPI  (Near  Kralak,  low)  _ UPPER  MISSISSIPPI  (Nw  KwAak.  lam)  _ UPPER  MISSISSIPPI  (New  Keokuk,  low) 
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Figure  5.17  Model-generated  monthly  precipitation,  runoff,  and  averaged  snow  mass  during  the  snow  season  for  the  present  climate 
[M(l)]  and  the  doubled  CO2  climate  [M(2)]  and  observed  [O]  precipitation  from  Legates  and  Willmolt  [1990]  and  runoff  from  the  US 
Geological  Survey  [1987]  for  the  Mississippi  River  at  Keokuk,  Iowa 


CENTRAL  MISSISSIPPI  (Near  St  Louis,  Mo.)  CENTRAL  MISSISSIPPI  (Near  St  Unis,  Mo.)  CENTRAL  MISSISSIPPI  (Near  St  Louie,  Mo.) 

- - — .  . . — — ~1  1 50 1 - - - -  ""  . '' 

Precipitation  Runoff  Sno»  1  x  CO, - 

2  X  CO,  •*.. 
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Figure  5  18  Model-generated  monthly  precipitation,  runoff,  and  averaged  snow  mass  during  the  snow  season  for  the  present  climate 
[M(1)J  and  the  doubled  C02  climate  [M(2)]  and  ouserved  [O]  precipitation  from  Legates  and  Willmoit  [1990]  and  runoff  from  the  US 
Geological  Survey  [1987]  for  the  Mississippi  River  at  St  Louis,  Missouri  and  Vicksburg,  Mississippi. 
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Geological  Survey  [1987]  for  the  Illinois  and  Ohio  rivers 
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Chapter  6.  Summary  and  Conclusions 

The  purpose  of  this  thesis  was  to  obtain  a  better  understanding  of  annual  and 
monthly  precipitation  and  river  runoff  as  it  is  modeled  in  several  simulations  for  the 
present  and  doubled  CO2  climates  using  the  GISS  GCM  of  Hansen  et  al.  [1983]  The 
river  basins  are  divided  into  categories  depending  on  precipitation  rates  [Kuhl  and  Miller , 
1992],  Observed  precipitation  was  compared  using  two  studies.  Shea  [1986]  and  Legates 
and  Willmott  [1990],  For  the  majority  of  the  river  basins,  annual  precipitation  was  higher 
(10%)  in  Legates  and  Willmott  [1990],  however,  the  monthly  variation  was  very  similar 
between  the  two  studies.  Therefore,  observed  precipitation  was  taken  from  Legates  and 
Willmott  [1990],  which  contained  more  observations  and  was  the  more  recent  of  the  two 
studies. 

The  mean  annual  and  total  basinwide  monthly  runoff  studies  of  Russell  and  Miller 
[1990],  Miller  and  Russell  [1992],  and  Kuhl  and  Miller  [1992]  were  extended  by 
allowing  river  runoff  to  be  routed  between  grid  boxes  within  a  river's  basin  using  the 
routing  scheme  of  Miller  et  al.  [1992],  The  timing  of  the  model -generated  total  river 
runoff,  Rtot,  without  the  routing  scheme  is  primarily  based  on  the  model's  precipitation 
This  was  clearly  illustrated  in  the  Congo  River  in  which  precipitation  and  runoff  maxima 
occurred  in  the  same  month.  However,  with  the  routing  scheme  model-generated  runoff 
at  the  river  mouth,  Rm,  was  lower  in  the  maximum  precipitation  month  and  higher  in  the 
following  months,  indicating  that  water  that  contributed  to  the  peak  in  Rtot  reached  the 
mouth  in  the  following  months  similar  to  observed  river  runoff. 

A  common  characteristic  of  the  observed  runoff  in  rivers  at  high  latitudes  in  the 
northern  hemisphere  is  a  peak  occurring  in  spring  and  early  summer  This  peak  is  due  to 
the  melting  of  snow  and  ice  The  new  routine  scheme  allows  water  to  move  downstream 
and  simulate  the  magnitude  and  timing  of  the  runoff  peak  at  the  river  mouth  In  this 
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particular  formulation  of  the  routing  scheme,  Q  and  Ay  of  Eqn  2.10  are  constants  The 
flow  rate  does  depend  on  topography  and  the  distance  across  the  grid  box,  As. 

As  noted  in  Miller  el  al.  [1992],  H  is  a  function  of  grid  resolution  It  is  also  likely 
to  be  dependent  on  other  characteristics  of  the  grid  box.  An  alternative  formulation  for 
Eqn.  2.10  given  by  Miller  et  al.  [1992]  can  eliminate  or  reduce  the  dependence  on  grid 
resolution.  However,  the  effects  of  the  physical  characteristics  of  a  particular  grid  box 
could  serve  as  a  basis  for  parameterizing  fi  within  each  grid  box  so  that  it  would  no 
longer  be  a  globally  uniform  constant  Ultimately  the  objective  is  to  accurately  route  river 
runoff  through  the  basin  with  routing  coefficients  based  on  the  simplest  possible 
parameterizations  of  Q 

The  model-generated  runoff  in  the  dry  regions  is  too  large.  The  combination  of 
too  little  evaporation  or  percolation  into  the  soil  and  too  much  model-generated 
precipitation  contribute  to  the  excess  computed  runoff  in  the  dry  regions  Also  as 
discussed  in  Miller  and  Russell  [1992],  the  model  does  not  allow  runoff  to  evaporate  or 
percolate  into  the  soil  as  it  moves  downstream  from  one  grid  box  to  another 

River  runoff  within  the  Mississippi  river  basin  was  examined  Although  the  model¬ 
generated  annual  precipitation  was  within  5%  of  the  observed  precipitation,  monthly 
precipitation  was  not  as  close  This  affected  the  model-generated  monthly  river  runoff 
Model-generated  precipitation  was  too  high  for  the  Missouri  River  and  too  low  for  the 
Ohio  River.  There  also  seemed  to  be  too  much  model-generated  snow  in  the  Missouri 
River  basin  and  not  enough  for  the  Ohio  River  basin  This  problem  can  be  caused,  in  part, 
by  the  resolution  of  the  grid  box.  Precipitation  and  snow  were  interpolated  from  a  4°  X  5° 
grid  box  in  which  precipitation,  snow,  and  runoff  were  assumed  equal  for  the  entire  grid 
box  In  the  case  of  the  Missouri  River,  the  western  edge  of  the  basin  borders  the  rocky 
mountains  with  only  half  of  the  4°  X  5°  grid  box  covering  the  Missouri  River  basin  It  is 
possible  that  most  of  the  precipitation  and  snow  accumulates  in  the  western  part  of  the  4° 
X  5°  grid  causing  higher  amounts  of  precipitation  and  snow  to  occur  in  the  Missouri  river 
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basin  after  interpolation  Problems  in  resolutions  near  mountainous  regions  can  cause 
erroneous  amounts  of  precipitation  and  runoff  to  occur  in  a  river  basin  Other  problems 
which  occur  in  mountainous  regions  are  grid  boxes  that  contain  more  than  one  basin 
Water  from  the  entire  grid  box  is  assumed  to  flow  in  one  direction,  thus,  some  portion  of 
the  river  runoff  may  enter  the  wrong  river  basin  Finer  resolution  models  could  be  used  to 
avoid  this  problem 

Tnis  study  extended  the  study  of  Miller  and  Russell  [1992]  to  examine  seasonal 
variations  in  river  runoff  in  a  doubled  C02  climate  They  found  that  for  the  doubled  C02 
climate,  mean  annual  runoff  increases  for  27  of  the  30  rivers  examined  The  annual  runoff 
increased  in  all  the  high  latitudes,  with  increases  averaging  approximately  28%  This  is 
consistent  with  other  studies  which  show  increasing  runoff  at  high  latitudes  for  a  doubled 
C02  climate  [Manabe  and  Stauffer,  1980,  Rind,  1988,  Mitchell,  1989,  Stauffer  et  ai. 
1989] 

Precipitation  and  temperatures  from  GCMs  have  been  used  in  hydrologic  models 
(off-line  method)  to  predict  regional  runoff  changes  that  would  accompany  global 
warming  [Gleick,  1987,  Flaschka,  et  ai,  1987,  Rons,  1989,  Lettenniaier  and  Gan,  1990. 
and  Ayers  et  al ,  1990]  Gleick  [1987]  and  Lettenmaier  and  Gan  [1990]  found  that  the 
principal  physical  mechanism  affecting  increased  winter  runoff  and  decreased  summer 
runoff  was  a  decrease  in  snow  as  a  proportion  of  the  winter  precipitation,  similar  to  the 
results  found  in  this  study 

The  most  interesting  changes  in  seasonal  runoff  occurred  in  the  higher  latitudes 
where  snow  melt  is  an  important  component  of  river  runoff  The  timing  of  the  snow  melt 
depends  on  temperature  Hence,  the  delays  for  the  present  climate  could  be  related  to 
temperature  as  well  as  the  routing  scheme  This  will  affect  the  timing  in  a  doubled  C02 
climate  as  well  Precipitation  increases  and  the  spring  runoff  peak  occurs  earlier  for  all 
high  latitude  rivers  Of  the  nine  high  latitude  basins  examined,  five  river  basins  in 
northeastern  Asia  (the  Yenesei,  Lena,  Kolyma,  Indigirka,  and  Amur)  also  show'  an 
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increase  in  the  magnitude  of  the  runoff  peak  This  region  is  the  only  region  to  show  an 
increase  in  snow  mass  in  a  doubled  C02  climate.  The  geographical  location  and 
climatology  of  this  region  indicate  that  winter  temperatures  do  not  increase  there 
sufficiently  to  reduce  the  winter  snowfall 

The  model-generated  river  runoff  decreases  36%  in  a  doubled  C02  climate  for  the 
Niger  River,  despite  little  change  in  the  annual  precipitation  between  the  present  and 
doubled  C02  climates  This  must  be  due  to  increased  evapotranspiration  which  would 
increase  more  in  regions  of  high  temperatures  because  of  the  non  linearity  of  the  Clausius 
Clapeyron  equation. 

Although  changes  in  river  runoff  were  obtained  for  a  doubled  C02  climate,  further 
improvements  in  the  model  are  needed  to  increase  the  confidence  in  the  results.  The 
limited  three-year  and  five-year  simulations  should  be  extended  to  reduce  the  chance  of 
temporal  errors.  Finer  resolutions  in  defining  river  basins  are  needed,  especially  in  areas 
which  border  mountainous  regions.  Also,  grid  boxes  need  to  be  divided  so  that  the  river's 
drainage  area  is  equal  to  the  appropriate  percentage  of  a  grid  box's  runoff 

The  conclusions  about  river  runoff  in  the  present  and  doubled  C02  climates 
depend  on  the  model’s  ability  to  simulate  the  hydrologic  cycle,  which  depends  o.  the 
model's  parameterizations  of  land-atmosphere  interactions.  It  is  essential  for  hydrologists 
and  climate  modelers  to  develop  the  best  possible  parameterizations  of  land-atmosphere 
interactions  within  GCMs  The  confidence  in  the  model's  ability  to  generate  snow  for  the 
present  climate  is  hard  to  evaluate  because  observed  snow  was  not  obtained  for 
comparison  with  model  snow.  Rind  et  al.,  [1990]  conclude  that  drought  intensification  is 
understated  in  a  doubled  C02  climate  for  several  GCMs,  including  the  GISS  model  used 
for  this  study,  because  of  their  failure  to  show  extensive  soil  moisture  reductions  This  is 
due  to  unrealistic  simulations  of  the  land  surface  A  new  soil-moisture  storage  scheme  has 
been  developed  by  Abramopoulos  et  al.,  [1988]  to  replace  the  simplified  two-layer  storage 
scheme  used  here,  but  the  new  scheme  has  not  been  implemented  fully  into  the  GCM 
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This  scheme  also  includes  groundwater.  A  primary  concern  of  climatic  modelers  must  be 
to  obtain  accurate  precipitation,  evapotranspiration,  and  soil  moisture  storage.  River 
runoff  provides  a  useful  diagnostic  for  examining  parameterizations  of  these  processes 
The  ability  to  predict  changes  in  river  runoff  are  essential  for  forecasting  future  water 
resource  needs. 
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APPENDIX:  Tables  of  monthly  statistics  for  each  river  basin 


River  runoff  (Km3)  and  precipitation  (Km3)  for  the  following  rivers: 


Wet 

Moderatelv  Wet 

Dm 

High  Latitude  Drv 

Amazon 

Yangtze 

Indus 

Yenesei 

Congo 

Mississippi 

Tigris-Euphrates 

Lena 

Orinoco 

LaPlata  (Parana) 

Yellow 

Ob 

Mekong 

St  Lawrence 

Colorado 

Amui 

Magdalena 

Danube 

Murray 

Mackenzie 

Sao  Francisco 

Columbia 

Yukon 

Zambesi 

Sevemay  Dvina 

Fraser 

Kolyma 

Nile 

Indigirka 

Niger 

Otos  -  Observed 

l(tot)  -  Model-generated  runoff  averaged  over  the  entire  river  basin 

H(m)  -  Model-generated  runoff  at  the  mouth  of  the  river  basin 

1XC02  -  Present  climate 

2XC02  -  Doubled  C02  climate 

848F  -  Model  simulation 

C863  -  Model  simulation 

1186  -  Model  simulation 

A51M  -  Model  simulation 

9471  -  Model  simulation 


Observed  runoff  is  from  UNESCO  [1969,  1974,  1985],  observed  precipitation  is  from 
Legates  and  Willmotl  [1990]  and  Shea  [1986],  and  areas  are  from  Millmatt  and  Meade 
[1983]  and  UNESCO  [1969,  1974,  1985], 
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